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FOREWORD 


This  summary  report  describes  research  performed  on  a  program 
sponsored  by  the  Office  of  Naval  Research,  Department  of  the  Navy, 
under  contract  N00014-74-C-0241 .  The  research  was  performed  under 
the  general  technical  direction  of  Dr.  Robert  Pohanka  and 
Dr.  Arthur  M.  Diness,  Metallurgy  Branch,  Office  of  Naval  Research. 

The  report  covers  work  performed  during  the  period  of 
February  1,  1980  to  August  31,  1982.  The  writers  are  pleased  to 
acknowledge  the  contributions  of  their  associates  at  Ceramic 
Finishing  Company. 


I.  INTRODUCTION 


Research  on  this  contract  which  extended  from  1974  to  1982  has 
included  the  following  topics  related  to  localized  impact  damage: 

1.  Development  of  a  comprehensive  analytical  model  for 
elastic-plastic  impact  by  particles  for  calculating 
the  contact  radius,  indentation  radius,  impact  force, 
penetration  time,  depth  of  damage,  remaining  strength, 
contact  time  and  coefficient  of  restitution. 

2.  Localized  impact  damage  in  stressed  members  consisting 
of  glass,  alumina  and  silicon  nitride. 

3.  Branching  of  Hertzian  cracks. 

4.  Energy  loss  mechanisms  during  impact. 

5.  Impact  on  a  viscous  medium  (glass). 

6.  Impact  damage  in  transformation  toughened  zirconia. 

7.  Comparison  of  static  and  impact  damage  in  ZnS. 

8.  Load-radial  crack  size  relations  in  ZnS. 

9.  Comparisons  of  contact  damage  in  hot  pressed  and 
CVD  ZnS. 

10.  Elastic  recovery  at  static  loading  and  impact 
indentations  in  MgF2* 

11.  Effect  of  mechanical  impedance  on  the  distribution 
of  stress  wave  energy  during  impact. 


12.  Methods  of  improving  the  impact  damage  resistance  of 
ceramics  including  low  elastic  midulus  surface  layers, 
friable  surface  layers  and  compressive  surface  layers 
(silicon  nitride,  silicon  carbide,  alumina,  stc.). 

In  addition  to  the  above  topics  which  directly  involve  contact 
damage,  the  following  supporting  research  topics  were  investigated: 

1.  Crack  branching  in  ceramics. 

2.  Strengthening  by  heat  treatments. 

3.  Strengthening  by  compressive  surface  stresses  including 
preparation  of  a  monograph  on  this  topic. 

4.  Comparison?;  of  flaws  at  fracture  origins  in  various 
ceramics . 

5.  Crack  growth  from  small  flaws  in  larger  grains  in  alumina. 

The  results  of  these  investigations  show  that,  at  least  in  some 

cases,  reasonable  predictions  of  crack  propagation  and  strength 
degradation  are  possible.  The  impact  damage  resistances  of  a  variety 
of  ceramics  can  be  improved  by  surface  treatment.  Increasing  the 
fracture  toughness  (K^)  of  the  underlying  material  is  also 
effective.  The  roles  of  hardness  and  Young's  modulus  are  more 
complex  in  the  impact  case  compared  with  the  static  case  because  of 
the  effect  of  these  properties  on  the  impact  load  and  stress 
wave  energy. 

The  present  report  consists  of  five  papers  which  are  presented 
as  separate  sections.  A  list  of  papers,  reports  and  presentations 


follows. 
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CONTACT  DAMAGE  IN  HOT  PRESSED  AND  CVD  ZINC  SULFIDE 


by 


H.  P.  Kirchner 
J.  A.  Tiracorda 
T.  J.  Larchuk 


ABSTRACT 


Hot  pressed  zinc  sulfide  was  indented  by  spherical  indenters 
under  static  and  impact  loading.  The  static  and  impact  damage 
characteristics  were  compared  with  earlier  results  for  CVD  zinc 
sulfide  showing  that  hot  pressed  zinc  sulfide  is  less  resistant  to 
contact  damage  than  CVD  zinc  sulfide. 
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Polycrystalline  zinc  sulfide  can  be  prepared  in  a  transparent  form 

suitable  for  use  in  optical  components  by  at  least  two  methods,  hot  pressing 

(HP  ZnS)  and  chemical  vapor  deposition  (CVD  ZnS) .  Evans  and  Wilshaw^ 

characterized  static  contact  damage  at  spherical  indentations  in  CVD  ZnS 

and  other  ceramics  and  determined  that  crack  extension  depends  primarily 

on  the  indentation  radius  and  the  ratio  of  hardness  to  toughness.  Shockey, 

(2) 

Dao,  and  Curran  characterized  the  impact  damage  at  spherical  indentations 

in  CVD  ZnS  and  correlated  the  damage  with  the  plastic  impression  strain. 

(3-5) 

Kirchner,  Richard  and  Larchuk  characterized  static  and  impact  damage 

at  spherical  indentations  in  CVD  ZnS  and  used  the  resulting  data  to  evaluate 
a  theory  of  elastic-plastic  impact,  to  correlate  the  radial  crack  length 
with  load  and  indentation  radius,  and  to  demonstrate  the  importance  of 
strain  rate  effects.  Van  der  Zwaag,  Hagan  and  Field^  investigated  deform¬ 
ation  mechanisms  at  static  contacts  in  CVD  ZnS.  Other  investigators^7  10^ 
have  studied  water  drop  Impact  on  ZnS.  In  these  investigations,  contact 
damage  in  HP  and  CVD  ZnS  has  not  been  compared. 

The  HP  ZnS* **  was  purchased  as  a  cylindrical  plate  3  in.  diam.  x  6mm 
thick.  This  plate  was  indented  on  the  large  surfaces  which  are  perpendicular 
to  the  hot  pressing  direction  by  static  and  impact  loading  using  spherical 
glass  indenters,  3mm  in  diameter  .  The  static  loads  ranging  from  111-890N 
were  applied  by  an  Instron  testing  machine.  The  HP  ZnS  plate  was  placed  on 
the  load  cell  and  the  spheres  were  fastened  to  the  crosshead.  A  compressed 
gas  (N2)  gun  was  used  to  accelerate  the  spheres  for  the  impact  tests. 

The  plate  was  mounted  perpendicularly  to  the  path  of  the  sphere.  The 
impact  and  rebound  velocities  were  measured  from  multiple  images  of 


*  IRTRAN  2,  Eastman  Kodak  Co.,  Rochester,  N.Y.  14650 

**  Glass  spheres,  No.  3000,  Walter  Stern,  Inc.,  Port  Washington,  N.Y. 


the  sphere  on  a  photograph  obtained  by  illuminating  the  sphere  at  known 
time  intervals  using  a  stroboscope  and  recording  the  images  using  a  camera 
with  an  open  shutter.  The  velocities  of  the  glass  spheres  ranged  from 
18-73  m-s-1. 

The  damage  was  characterized  by  optical  and  scanning  electron  micro¬ 
scopy  and  the  results  were  compared  with  earlier  results  obtained  for  CVD 
ZnS  indented  on  the  surface  perpendicular  to  the  deposition  direction. 

I 

Comparison  of  microstructure  and  properties  of  HP  and  CVD  Zns 

The  microstructures  of  HP  and  CVD  ZnS  differ  because  of  the  differences 
in  the  forming  method.  The  HP  ZnS  has  a  grain  size  of  1-5  pm  which  does 
not  vary  significantly  with  orientation  in  the  plate.  The  CVD  ZnS  is 
deposited  in  the  form  of  columnar  grains  oriented  roughly  parallel  to  the 
deposition  direction.  These  grains  have  an  aspect  ratio  of  8:1  and  are 
up  to  perhaps  50  pm  in  length.  In  the  plane  perpendicular  to  the  deposition 
direction  the  grain  size  is  quite  variable  ranging  from  about  1-20  pm  and 
the  grains  are  approximately  equidimensional. 

The  properties  of  the  HP  and  CVD  ZnS  are  compared  in  Table  I  using 
data  from  various  sources.  The  of  the  HP  ZnS  was  measured  by  the  single 

edge  notched  beam  method  for  a  crack  propagating  parallel  to  the  hot 

k 

press* ng  direction  yielding  0.44  MPam  ,  a  value  substantially  lower  than 
the  literature  values  for  CVD  ZnS.  This  result  suggests  that  cracks  will 
propagate  farther  in  HP  ZnS  than  in  CVD  ZnS.  This  crack  propagation 
direction  is  relevant  to  propagation  of  the  radial  and  median  cracks.  Com¬ 
parisons  of  the  diameters  of  static  indentations  in  the  HP  and  CVD  ZnS 


_2 

show  that  the  hardnesses  are  approximately  equal  (1.9  GNm  ). 


TABLE  I 

Comparison  of  Properties  of  HP  and  CVD  ZnS 
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HP  ZnS* ** *** **** 

ieic 

CVD  ZnS 

Microstructure 

Grain  size,  um 

1-5 

*** 

(c,d,e,see  note  ) 

2-5  in  growth  plane 

9-30  X  "  " 

Crystal  habit 

- 

columnar  grains 

Crystal  symmetry 

cubic  ^ 

cubic 

Critical  stressx intensity 
factor,  MPa  m2 

**** 

0.44  ±  0.08 

1  <f> 

0.75(g) 

0.67(c) 

Density,  g-cm 

4.09(b) 

4.08(e) 

Elastic  properties 

Young's  modulus,  MPa 

96,550(b) 

74,500(e) 

Shear  modulus,  MPa 

- 

41,000(f) 

Poisson's  ratio 

- 

fp) 

0.30v  in  growth  plane 

0.35  X  " 

Modulus  of  rupture 
at  25°C,  MPa 

97.2(b) 

110.3(e) 

*  IRTRAN  2,  Eastman  Kodak  Co.,  Rochester,  N.Y.  14650. 

**  RAYTRAN  ZnS,  Raytheon  Research  Division,  Waltham,  Mass.  02154. 

***  References  d  and  e  indicate  wider  ranges  than  those  listed. 

****  Average  of  two  values. 

(a)  Personal  communication  with  Eastman  Kodak  Co.  (Mr.  Fairbanks). 

(b)  Manufacturer's  literature,  Eastman  Kodak  Co. 

(c)  J.  M.  Wimmer,  "Mechanical  and  Physical  Properties  of  Chemically 
Vapor-Deposited  (CVD)  Zinc  Sulfide,"  Air  Force  Materials  Laboratory 
Report  AFML-TR- 79-401 3  (1979). 

(d)  D.  A.  Shockey,  K.  C.  Dao  and  D.  B.  Curran,  "Nucleation  and  Growth  of 
Cracks  in  CVD  ZnS  under  Particle  Impact,"  SRI  International  Annual 
Report  Part  II,  Contract  N00014-76-C-0657  (April,  1979). 

(e)  Manufacturer's  literature,  Raytheon  Co. 

(f)  A.  G.  Evans  and  T.  R.  Wilshaw,  "Quasi-Static  Solid  Particle  Damage  in 
Brittle  Solids  I.  Observations,  Analysis  and  Implications,"  Acta  Met. 
24,939-956  (1976). 

(g)  S.  van  der  Zwaag,  J.  T.  Hagan  and  J.  E.  Field,  "Studies  of  Contact 
Damage  in  Polycrystalline  Zinc  Sulfide,"  J.  Mater.  Sci.  15  (12) 
2965-2972  (1980). 
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Comparison  of  static  and  impact  loading  damage  in  HP  and  CVD  ZnS 

The  expected  load  (P)  dependence  of  the  radial  crack  length  (c)  in  the 

h  (4) 

elastic-plastic  case  is  described  by  P  *  dc 2  where  d  is  a  contact  dimension 
Several  mechanisms  can  give  rise  to  this  relation  including  line  contact 
loading  or  wedging  where  d  is  the  length  of  the  contact ,  the  existence 
of  an  irreversibly  deformed  zone  under  the  contact  where  d  is  the  depth  of 
this  zone^"*’*^  and  the  existence  of  residual  stresses  where  d  is  the 
length  subject  to  stress  .  In  addition,  Niihara,  Morena  and  Hasselman^^ 
have  proposed  that  a  P  <*  c2  relation  applies  for  Palmquist  cracks. 

Indentations  and  cracks  formed  by  static  loading  in  HP  and  CVD  ZnS 
are  compared  in  Figure  1 ,  which  shows  that ,  compared  at  approximately  equal 
indentation  radii,  the  cracks  are  slightly  longer  in  HP  ZnS  than  in  CVD  ZnS 
and  also  shows  the  presence  of  a  lateral  crack  (the  circular  shadow)  in 
the  HP  ZnS.  The  load  threshold  for  lateral  crack  formation  is  lower  in 
HP  ZnS  than  in  CVD  ZnS.  Loads  of  500-600N  are  required  to  cause  lateral 
cracking  in  CVD  ZnS. 

The  extent  of  radial  cracking  in  HP  and  CVD  ZnS  is  compared  in  Figure 

2.  The  plots  show  that  the  variations  of  crack  length  with  load  in  both 

w 

materials  are  equally  well  represented  by  P  *  d^c  where  di  is  the  indentation 
diameter  and  that,  when  compared  on  this  basis,  the  cracks  are  longer  in 
HP  ZnS  than  in  CVD  ZnS.  At  high  loads  the  radial  cracks  in  HP  ZnS  are 
about  30%  longer  than  those  in  CVD  ZnS. 

Impacts  using  the  glass  spheres  yielded  damage  that  was  similar  in 
magnitude  to  the  damage  induced  by  static  loading  as  shown  in  Figure  3. 

For  specimens  impacted  at  approximately  equal  velocities,  the  radial  cracks 
are  longer  in  HP  ZnS  than  in  CVD  ZnS  as  was  the  case  for  static  loading. 


Static  Load 


O) 

6 

I  6.0 


4.0 


0  2.0  4.0  6.0  8.0  10.0 

(Radial  Crack  Length)1'2  -IO“2m,/2 

Figure  2.  Static  load  divided  hv  indentation  diameter  vs  square  root 
of  maximum  radial  crack  length  (for  both  CVI)  X.nS  and  Hr  ZnS 
statically  loaded  with  i  mm  diameter  glass  spheres). 


t\  =  286  pm,  c  =  792  pm. 


HP  ZnS  (impact  velocity  =  27  m*s  , 
c^  =  276  pm) 


CVD  Zns  (impact  velocity  =  28.1 
c  ^  =  82  pm) 


-1 

m-  s  , 


276  pm,  c^.  =  600  pm, 


Figure  3.  Damage  induced  by  impact  loading,  39X  (HP  ZnS  and  CVD 
ZnS  plates  with  3  mm  glass  spheres). 
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The  HP  ZnS  specimen  in  Figure  3  also  shows  a  partial  lateral  crack.  The 
differences  in  the  radial  crack  lengths  were  not  consistent  throughout  the 
velocity  range.  In  some  cases  the  radial  crack  lengths  in  HP  ZnS  were  much 
longer  than  would  normally  be  expected. 

The  circumferential  cracks  formed  in  the  impact  indentations  are  much 
better  defined  than  those  in  the  static  indentations. 

The  coefficients  of  restitution  (the  ratios  of  the  rebound  velocities 
to  the  impact  velocities)  were  calculated.  These  values  were  only  slightly 
different  for  the  two  materials  and  ranged  from  a  maximum  of  about  0.94 
at  the  low  end  of  the  velocity  range  to  a  minimum  of  about  0.76  at  the  high 
end  of  the  velocity  range.  The  variations  of  the  kinetic  energy  loss 
with  impact  velocity  were  compared.  Because  of  the  importance  of  plastic 
deformation  as  an  energy  absorption  mechanism,  differences  in  kinetic  energy 
loss  are  expected  to  indicate  differences  in  the  relative  dynamic  hardnesses 
of  the  materials.  The  slightly  higher  energy  losses  in  HP  ZnS  indicate 
that  the  dynamic  hardness  of  HP  ZnS  may  be  slightly  lower  that  that  of 
CVD  ZnS. 

The  log-log  plots  of  the  kinetic  energy  loss  vs.  impact  velocity  for 
both  HP  and  CVD  ZnS  are  essentially  linear  with  slopes  of  2.75  and  2.92 
respectively.  These  slopes  can  be  compared  with  a  value  of  2.6  predicted 
for  elastic  impacts  by  Hunter  This  result  provides  additional  evi¬ 

dence  that  the  presence  of  Impact  damage,  including  plastic  deformation  at 
indentations,  does  not  substantially  affect  these  slopes  except  at  transi¬ 
tions  where  new  energy  loss  mechanisms  become  Important. 

Therefore,  it  is  concluded  that  at  static  spherical  indentations  in 

HP  ZnS,  the  variation  of  the  radial  crack  length  with  load  is  well  repre¬ 
ss 

sented  by  P  «  d^  c  ,  The  region  in  which  this  relation  provides  a  good 
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fit  extends  to  c/r^  ratios  as  great  as  7.3,  much  larger  values  than  would 

Qg\ 

have  been  expected  based  on  previous  investigations  using  Vickers  indentersv  ' 
which  suggests  that  further  investigations  of  indenter  shape  and  stiffness 
might  be  fruitful.  At  high  static  loads  the  maximum  radial  crack  lengths 
in  HP  ZnS  are  about  30%  longer  and  the  numbers  of  radial  cracks  are  greater 
than  in  CVD  ZnS.  The  threshold  loads  for  lateral  crack  formation  are  much 
lower  in  HP  ZnS  than  in  CVD  ZnS. 

Under  comparable  impact  conditions  the  radial  crack  lengths  are 
greater  and  the  cracks  are  more  numerous  in  HP  ZnS  than  in  CVD  ZnS.  Also, 
crushing  is  more  extensive  in  HP  ZnS.  The  kinetic  energy  losses  during 
impact  are  only  slightly  greater  in  HP  ZnS  than  in  CVD  ZnS. 

Observations  of  longer  radial  cracks,  larger  numbers  of  radial  cracks, 
a  lower  threshold  load  for  lateral  crack  formation  and  increased  crushing 
in  HP  ZnS  indicate  that  HP  ZnS  is  more  susceptible  to  contact  damage  than 
CVD  ZnS.  This  difference  is  attributable,  at  least  in  part,  to  the  dif¬ 
ference  in  the  KjC  values. 
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ABSTRACT 

Contact  damage  was  induced  in  CVD  zinc  sulfide  by  static  and  impact 
loading  using  glass  spheres.  Calculations  based  on  the  relationship  be¬ 
tween  the  residual  depth  of  the  indentation  and  the  ratio  of  the  hardness 
to  Young's  modulus  were  used  to  estimate  the  dynamic  hardness  and  impact 
load.  Comparisons  of  applied  loads,  threshold  loads  for  crack  formation 
and  crack  lengths  for  static  and  impact  indentations  at  equal  indentation 
radii  show  that  predictions  of  impact  damage  using  data  from  static  load¬ 
ing  are  inadequate. 
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I.  INTRODUCTION 

A  common  approach  to  investigation  of  localized  impact  damage  has 
been  to  characterize  the  damage  at  static  indentations  induced  by  various 
loads  and  to  use  the  results  of  these  observations  to  predict  the  type 
and  extent  of  impact  damage,  assuming  that,  at  equal  loads,  identical  dam¬ 
age  is  induced  irrespective  of  the  method  of  loading.^  ^  Evaluation  of 
this  assumption  is  handicapped  by  the  difficulty  in  measuring  the  impact 
load.  However,  the  validity  of  the  assumption  can  be  evaluated  indirectly 
by  comparing  damage  at  static  and  impact  indentations  that  have  an  equal 
measurable  characteristic  such  as  the  indentation  radius.  Presumably,  if 
the  response  of  the  material  depended  only  on  the  maximum  load  and  not  on 
the  loading  rate,  each  characteristic  of  the  damage  observed  at  equal  in¬ 
dentation  radii  will  be  the  same  for  both  static  and  impact  loads.  The 
extent  of  the  observed  differences  will  indicate  the  degree  to  which  the 
response  differs  from  the  assumption. 

The  details  of  the  damage  induced  during  impact  are  important  because 
of  the  differing  roles  of  the  various  types  of  cracks  in  the  performance 
of  the  material.  For  example,  radial  cracks  are  considered  to  control  the 
strength  degradation.  Lateral  cracks  control  the  degradation  of  optical 
transmission  and,  combined  with  crushing,  control  the  erosion  rate. 
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II.  PROCEDURES 

When  a  sufficiently  high  load  is  applied  to  a  sphere  pressed  against 
a  flat  plate,  the  yield  point  of  the  plate  material  will  be  exceeded  and 
a  permanent  indentation  is  induced  (Figure  1).  As  the  load  is  removed, 
some  of  the  elastic  deformation  is  recovered  while  the  remainder  is  con¬ 
strained  by  the  irreversibly  deformed  material  resulting  in  residual 
stresses . 

In  this  investigation,  zinc  sulfide  plates  (25x25x7mm)  were  indented 
by  glass  spheres,  3mm  diameter,  under  various  loads  using  an  Instron  test¬ 
ing  machine.  The  materials  are  identified  and  the  properties  of  the 
plates  and  spheres  are  given  in  Table  I.  The  zinc  sulfide  consists  of 

(4) 

columnar  grains  with  a  high  degree  of  preferred  orientation  .  The  ob¬ 
served  grain  size  ranges  from  one  to  40  pm  depending  especially  on  the 
orientation  relative  to  the  deposition  direction.  The  indentation  radius 
(r^)  was  measured  in  the  plane  of  the  surface  as  indicated  in  Figure  1. 

The  indentation  radii,  numbers  of  radial  cracks  and  maximum  radial  and 
lateral  crack  lengths  were  measured,  using  a  calibrated  eyepiece  in  an 
optical  microscope  for  the  load  range  of  84  to  1059  N. 

Similar  plates  were  impacted  at  velocities  ranging  from  19  to  65  m-s 
by  glass  spheres  accelerated  by  a  gas  gun  using  compressed  N?  gas.  The 
impact  and  rebounding  velocities  were  measured  using  photographs  wi th 
multiple  stroboscopic  images  of  the  glass  spheres.  The  spheres  were 
coated  with  a  thin  layer  of  aluminum  to  improve  their  reflectivity.  The 
impact  damage  was  characterized  as  described  above  for  the  static  indenta¬ 
tions.  In  some  cases  the  profiles  of  the  static  and  impact  indentations 
were  compared,  using  a  profilometer.  The  results  were  analyzed  and  are 
presented  in  the  following  section. 


maximum  depth  (Z 
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*  Glass  spheres,  No.  3000,  Walter  Stern,  Inc.,  Port  Washington,  N.Y. 

**  Chemical  vapor  deposited  zinc  sulfide  plates,  Raytheon  Co., 
Waltham,  Mass. 
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III.  RESULTS  AND  DISCUSSION 

The  response  of  the  ZnS  plates  to  static  and  impact  loading  was 
elastic-plastic,  yielding  permanent  indentations  with  the  indentation 
radii  increasing  with  the  maximum  loads.  Ring  cracks,  radial  cracks  and 
lateral  cracks  were  observed  for  both  static  and  impact  loading  as  shown 
in  Figure  2  where  the  lateral  cracks  appear  as  the  shaded  areas  between 
the  radial  cracks.  When  the  damage  was  observed  in  cross  section,  a 
clear  zone  was  observed  beneath  the  contact.  This  zone  is  relatively 
transparent  because  of  the  absence  of  large-scale  cracks.  However,  this 

(4) 

zone  contains  dispersed  spots  that  have  been  identified  as  pores  . 

At  static  indentations,  the  load  thresholds  w«  re  84  N  for  radial 
crack  formation  and  899  N  for  lateral  crack  formation.  In  the  impact  tests 
measurable  indentations,  radial  cracking,  and  evidence  of  lateral  crack¬ 
ing  were  observed  at  the  lowest  velocity  tested  and  all  higher  velocities. 

Damage  at  static  and  impact  sites  was  compared  at  approximately 
equal  indentation  radii  (Figure  2).  When  the  indentation  radii  are  rela¬ 
tively  small,  as  in  Figures  2A  (static)  and  2B  (impact),  the  radial  cracks 
formed  at  the  impact  sites  are  fewer  in  number  and  smaller  than  those 
formed  at  the  static  sites.  However,  in  the  impact  case,  lateral  cracks 
initiate  at  smaller  indentation  radii  as  indicated  by  the  lighter  area 
to  the  left  of  the  indentation.  Also,  the  ring  cracking  is  much  more 
visible  at  the  impact  site.  At  larger  indentation  radii,  as  shown  in 
Figure  2C  (static)  and  2D  (impact),  both  the  radial  cracks  and  the  lateral 
cracks  are  larger  at  the  impact  site  than  at  the  comparable  static  site. 

The  lateral  cracking  at  the  impact  site  has  produced  chipping  to  the  right 


29 


of  the  indentation.  Also,  the  ring  cracking  within  the  impact  indentations 
has  become  much  more  pronounced,  causing  the  contours  of  these  indentations 
to  appear  somewhat  irregular,  compared  with  the  static  indentations. 

The  irregularities  observed  in  the  impact  indentations  were  investi¬ 
gated  by  profilometer  measurements.  Figure  3  shows  the  diametral  profiles 
of  static  and  impact  indentations,  compared  at  approximately  equal  inden¬ 
tation  radii.  It  should  be  noted  that  the  vertical  scale  is  ten  times  the 
horizontal  scale.  The  distinguishing  features  are: 

1.  The  impact  indentation  is  much  shallower  than  the 
static  indentation. 

2.  The  center  of  the  impact  indentation  is  separated  from 
the  outer  rim  by  substantial  circumferential  (ring) 
cracks,  one  of  which  is  visible  at  the  boundary  of  the 
dark  central  area  in  Figure  2D. 

3.  ihe  central  area  of  the  impact  indentation  has  been 
pushed  down,  apparently  as  a  result  of  shear  failure 
at  the  circumferential  cracks. 

4.  The  outer  portions  of  the  impact  profile  form  a  rim 
with  a  much  lower  slope  than  that  of  a  comparable  static 
indentation. 

These  features  indicate  that,  mechanically,  the  formation  of  impact 

indentations  is  much  more  complex  than  the  formation  of  static  indentations. 

(9) 

The  above  results,  together  with  recent  equations  of  Lawn  and  Howes  , 
were  used  to  estimate  the  dynamic  hardness  and  the  maximum  impact  load 
for  this  impact.  Lawn  and  Howes  expressed  the  elastic  recovery  for  elastic- 
plastic  indentations  by  Vickers  indenters  in  terms  of  the  hardness  to 
Young's  modulus  ratio  (H/E)  obtaining  the  relation 


Static  site  Impact  site 

q  3  400  ftm  s  396  yxm 

Load  *  936N  V0  *  61.6  m/s 

Zr  ■  30.4 /im  Zr  *  2l.6ft.m 


Scales: 


50^i.m 


500/i.m 


Note  difference 
between  horizontal 
and  vertical 
scales. 


Figure  3.  Diametral  profilometer  traces  of  static  and  impact 
indentations  (ZnS,  3  mm  diameter  glass  spheres). 
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where  Z  is  the  maximum  penetration  under  load,  Z  is  the  depth  of  the 
indentation  after  unloading,  H  is  the  hardness,  E  is  Young's  modulus, 
v  is  Poisson's  ratio,  4*  is  the  indenter  angle  and  and  y„  are  geometrical 
factors.  In  the  case  of  Vickers  indenters  'P  is  fixed  by  the  pyramidal 
shape  and  the  rigidity  of  the  diamond.  The  above  equation  might  be  ex¬ 
pected  to  apply  to  glass  spheres  as  indenters,  despite  the  different 
shape  and  the  elastic  deformation  if  an  appropriate  value  of  tan  V  can  be 

k 

estimated  in  each  particular  case.  Also,  Z  cannot  be  estimated  directly 
from  the  spherical  geometry  of  the  indenter  and  the  indentation  radius 
(r^)  because  of  the  elastic  deformation.  However,  7,  for  the  impact  case 
in  Figure  3  can  be  estimated  from  measured  values  of  the  coefficient  of 

(9) 

restitution  (e)  using 


(5)  * 

In  this  case  e  =  0.79  and  Z  =  40.75  pm.  Then,  noting  that  r^  is  approxi¬ 
mately  the  same  in  both  the  static  and  impact  cases,  it  can  be  assumed  that 
* 

Z  is  approximately  the  same  in  both  cases,  and  that  =  yh  so  that  tan  V  can 
be  estimated  from  equation  (1)  for  the  static  case  yielding  Y  =  83.8° 
where  Y  in  this  case  is  an  apparent  indenter  angle.  Then,  assuming  that 
tan  Y  is  approximately  the  same  in  both  the  static  and  impact  cases  the 


dynamic  hardness  can  be  estimated  from  equation  (1)  for  the  impact  case 
yielding  3.1  GPa.  This  dynamic  hardness  implies  an  impact  load  of 
1500  N  which  can  be  compared  with  the  static  load  of  936  N  at  approxi¬ 
mately  the  same  indentation  radius  and  an  impact  load  of  1800  N  estimated 
previously using  the  Hertzian  equations  that  assume  elastic  conditions 
Because  the  plastic  deformation  reduces  the  maximum  load  one  would  expect 
the  load  in  the  elastic-plastic  case  to  be  less  than  that  in  the  elastic 
case  as  is  indicated  above. 

One  can  speculate  that  the  boundary  formed  by  the  wide  circumfer¬ 
ential  crack  represents  a  boundary  between  a  central  region  in  which  force 
tending  to  cause  radial  flow  of  the  ZnS  along  the  interface  between  the 
ZnS  and  glass  are  insufficient  to  overcome  the  frictional  forces and 
a  circumferential  band  in  which  the  frictional  forces  are  overcome  so  that 
there  is  movement  of  ZnS  radially  along  this  interface. 

The  extent  of  damage  under  the  surface  was  observed  in  detail  by 
viewing  the  damage  through  the  polished  sides  of  the  plates  as  shown  in 
Figure  4  in  which  the  radial  cracks  are  evident.  The  sketches  in  Figun  5 
present  details  of  cracks  observed  in  this  manner.  Out-of-plane  radial 
cracks  were  omitted  for  clarity.  The  pairs  of  static  and  impact  specimens 
have  comparable  indentation  radii  and  are  arranged  in  order  of  increasing 
indentation  radius.  The  cross-section  in  Figure  5A  represents  damage 
at  a  small  static  indentation.  The  damage  consists  of  radial  (Palmquist) 
cracks  extending  from  the  edges  of  the  indentation  and  a  clear,  ir¬ 
reversibly  deformed  zone  under  the  indentation.  The  radial  cracks  to  not 
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Figure  4 
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Side  view  of  damage  under  the  surface  at  a  static 
indentation  (load  =  708  N,  r ^  =  348  pm). 
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Static  site 


Static  site 


Lood  *  296N 
r(  *  247  pm 


Load  ■  708N 
fj  *  348  pm 


A 
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Impact  site 


V0  s  19  7  m/s 
f(  1  24  3  pm 


V0  «  40  8  m/s 
r  :  351  p  m 


B 


0 


Figure  5.  Comparisons  of  cracks  formed  ac  static  and  Impact  sites 
with  comparable  indentation  radii  (ZnS,  3  mm  diameter 
glass  spheres). 
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extend  downward  into  the  plate  appreciably  until  they  are  beyond  the 

irreversibly  deformed  zone.  The  radial  cracks  are  not  necessarily  coplanar 

and  they  do  not  bridge  the  gap  under  the  indentation. 

(2) 

Evans  and  Wilshaw  observed  that  radial  cracks  form  on  loading  in 

ZnS.  This  observation  is  confirmed  by  observation  of  radial  cracks  that 

extend  into  the  indentation  (Figure  2A) .  If  the  radial  cracks  formed  on 

unloading,  as  a  response  to  the  residual  stresses,  the  radial  cracks  would 

originate  at  the  boundary  of  the  irreversibly  deformed  zone. 

The  Boussinesq  and  Michell  stress  distributions  for  a  point  load  and 

the  Hertzian  stress  distribution  for  a  blunt  load  at  the  surface  show  that, 

in  the  absence  of  plastic  deformation,  the  out-of-piane  hoop  stresses  are 

compressive However,  Perrott  has  shown  that  when  the  plastic 

zone  radius  is  1.65  times  the  indentation  radius  the  out-of-plane  hoop 

stresses  at  the  boundary  of  an  indentation  become  tensile  at  the  surface. 

Also,  for  an  extreme  case  in  which  a  rigid  sphere  is  pressed  to  the 

depth  of  its  radius  into  zinc  sulfide,  numerical  analysis  of  the  elastic 

(2) 

stress  field  by  Evans  and  Wilshaw  has  shown  that  the  out-of-plane  hoop 
stresses  are  tensile.  Therefore,  the  observation  that  radial  cracks  form 
on  loading  in  zinc  sulfide  can  be  attributed  to  the  extensive  plastic 
deformation  of  this  relatively  soft  material. 

The  impact  site  (Figure  5B)  differs  from  the  comparable  static  site 
in  several  ways.  The  impact  site  has  only  one  radial  crack,  which  is 
smaller  than  the  largest  radial  crack  at  the  static  site.  Also,  lateral 
cracking  is  present  and  the  irreversibly  deformed  zone  is  smaller.  These 
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observations  are  consistent  with  the  earlier  observation  that  the  impact 
indentations  are  shallower  than  static  indentations  at  comparable  inden¬ 
tation  radii.  Therefore,  one  would  expect  lower  tensile  out-of-plane  hoop 
stress  during  loading  and  lower  residual  stresses  during  unloading  which 
would  result  in  fewer  and  smaller  radial  cracks.  Also,  because  the  dynamic 
hardness  is  greater  than  the  static  hardness  in  ZnS,  the  presence  of 

lateral  cracking  at  the  impact  site  is  consistent  with  the  observation  of 

(2) 

Evans  and  Wilshaw  that  lateral  cracking  is  more  extensive  in  harder 
materials . 

At  the  next  larger  indentations  (Figures  5C,D),  Lhe  radial  cracks  at 
the  static  indentation  are  much  larger  than  those  in  Figure  5A,  and  they 
have  extended  so  that  they  almost  meet  under  the  indentation.  Otherwise, 
the  damage  is  similar  to  that  in  Figure  2A.  The  impact  damage  shows 
larger  radial  and  lateral  cracks,  as  expected.  In  addition,  under  the 
clear  zone,  there  are  cracks  extending  deeper  into  the  ZnS  perpendicular 
to  the  surface  that  have  been  labelled  as  median  cracks.  These  cracks  are 
not  coplanar  and  the  vertical  line  extending  along  the  impact  axis  repre¬ 
sents  the  intersection  of  the  cracks.  Median  cracks  are  frequently  observed 
at  static  indentations  under  sharp  indenters.  As  the  load  is  removed, 
radial  cracks  may  extend  in  response  to  residual  stresses  and  merge  with 
existing  median  cracks.  Apparently,  this  has  occurred  in  the  present  case. 

At  the  largest  indentations  (Figure  5E,F),  the  radial  cracks  at  the 
static  indentation  are  even  larger  and  they  have  extended  under  the  inden¬ 
tation  to  completely  close  the  gap  to  form  an  approximately  semicircular 


crack.  A  large  lateral  crack  was  observed  close  to  the  surface.  The  im¬ 
pact  damage  shows  larger  non-coplanar  radial  cracks  and  larger  lateral 
cracks  compared  with  those  in  Figure  5D.  Again,  median  cracks  have  merged 
with  radial  cracks.  Lateral  cracks  are  observed  at  several  levels,  all 
farther  from  the  surface  than  the  lateral  crack  at  the  comparable  static 
site. 

The  numbers  and  lengths  of  radial  cracks  in  static  and  impact  loaded 
specimens  are  compared  in  Figures  6  and  7.  Radial  cracks  at  static  inden¬ 
tations  are  observed  at  smaller  indentation  radii  and  larger  numbers  of 
radial  cracks  are  formed  at  a  given  indentation  radius  than  at  impact 
indentations  (Figure  6).  Apparently,  the  more  extensive  irreversible 
deformation  at  the  static  indentations  changes  the  sign  of  the  out-of¬ 
plane  hoop  stresses  from  compressive  to  tensile  at  smaller  indentation 
radii  than  in  the  impact  case  leading  to  greater  tensile  hoop  deformation 
and  the  greater  number  of  cracks.  Also,  the  extensive  irreversible  de¬ 
formation  may  induce  substantial  residual  stresses. 

Radial  cracks  at  static  indentations  are  longer  than  those  at  impact 
sites  for  indentation  radii  up  to  about  400  urn  (Figure  7).  At  greater 
radii  the  radiaL  cracks  at  impact  sites  are  longer.  Apparently,  at  small 
indentation  radii,  the  deeper  indentations  and  greater  irreversible  de¬ 
formation  observed  at  the  static  sites  changes  the  sign  of  the  out-of¬ 
plane  hoop  stresses  at  lower  loads  and  results  in  greater  residual  stresse 
leading  to  longer  radial  cracks  than  those  at  comparable  impact  indenta¬ 
tions.  However,  at  large  indentation  radii,  the  larger  loads  resulting 
from  increased  dynamic  hardness  during  impact  become  dominant  and  the 
cracks  become  larger  than  those  at  comparable  static  indentations. 
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Figure  7.  (’rack  lengths  for  static  and  impact  radial  and  lateral 
cracks  vs.  indentation  radius  (ZnS,  3  mm  diameter 
glass  spheres). 
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IV.  SUMMARY  AND  CONCLUSIONS 

ZnS  plates  were  statically  indented  under  various  loads  and  impacted 
at  various  velocities  using  3mm  diameter  glass  spheres.  The  contact 
damage  was  characterized  and  the  damage  at  the  static  and  impact  sites 
was  compared  yielding  the  results  described  in  the  previous  section.  The 
results  show  that  when  compared  at  equal  indentation  radii,  there  are 
substantial  differences  between  static  and  impact  damage. 

Damage  at  the  static  and  impact  indentations  included  the  irreversi¬ 
bly  deformed  zone  directly  under  the  indentation  and  ring,  radial,  and 
lateral  cracks.  The  scale  of  the  damage  at  the  static  indentations 
increased  with  increasing  loads  on  the  glass  spheres.  An  overall  compari¬ 
son  of  the  static  and  impact  damage  at  approximately  equal  indentation 
radii  shows  that  the  scale  of  the  overall  impact  damage  is  approximately 
the  same  as  that  of  the  static  indentation  damage.  If  one  assumes  that 
the  sizes  of  the  cracks  mainly  depend  on  the  maximum  loads  applied  at 
the  contact,  this  observation  appears  to  justify  the  assumption  of  equal 
loads  at  equal  indentation  radii  for  static  and  impact  indentations. 
However,  the  calculations  and  microscopic  observations  show  that,  at 
equal  indentation  radii  the  impact  loads  are  substantially  higher  than 
the  static  loads  and  there  are  substantial  differences  in  the  indentation 
radius  thresholds  for  crack  formation  and  the  numbers  and  sizes  of  the 
various  types  of  cracks  at  static  and  impact  indentations.  These  differ¬ 
ences  are  important  in  the  performance  of  the  material.  For  example,  pre¬ 
dictions  based  on  the  static  damage  may  overestimate  the  extent  of  radial 
cracking  and  strength  degradation  at  low  impact  velocities.  Also,  they 
may  underestimate  the  extent  of  lateral  cracking  and,  therefore,  under¬ 
estimate  the  degradation  of  optical  transmission  and  erosion  rate. 
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ABSTRACT 

HP  ZnS  was  indented  by  spherical  indenters  of  various  radii.  The 

relationship  of  the  load  (P)  and  radial  crack  length  (cr)  was  better 

u 

modeled  by  P  «  d.  c  2  where  d.  is  the  indentation  diameter  than  hv 
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I.  INTRODUCTION 


(1  2)  (2-6) 

Elasticv  ’  and  elastic-plastic'  '  analyses  have  been  used  to 
relate  the  extent  of  crack  propagation  to  the  applied  load  at  contacts 
with  various  types  of  indenters.  The  resulting  relations  can  be  used 
to  estimate  the  extent  of  crack  propagation,  strength  degradation  and 
erosion  resulting  from  damage  induced  during  static  indentation^ 
localized  impact  and  grinding^ ^  The  available  relations 

differ  mainly  in  the  extent  to  which  they  take  account  of  the  character¬ 
istics  of  the  particular  contact,  including  the  effect  of  residua] 
stresses.  In  the  region  far  removed  from  the  contact  (the  far  field), 

the  effects  of  the  particular  contact  characteristics  are  negligible 

3/  , 

and  relations  of  the  form  load  P  «  c  '  where  c  is  the  crack  dimension 
can  be  used.  However,  in  a  wide  range  of  cases,  the  effects  of  the 
contact  characteristics  are  important  and  must  be  accounted  for.  The 
contact  characteristics  may  be  important  when  the  contact  area  is 
large,  when  there  is  extensive  irreversible  deformation  at  the  contact, 
in  the  presence  of  residual  stresses,  and  in  cases  of  line  contact 
loading  or  wedging.  In  this  investigation,  the  indenter  radii  and 
applied  loads  were  varied  and  the  resulting  observations  were  used 
to  evaluate  the  effect  of  contact  characteristics  associated  with 


indenter  radius  on  the  extent  of  crack  propagation  in  zinc  sulfide. 


The  original  elastic  relation  for  the  far  field  was  derived  by 

(1  2) 

Roesler  ’  for  the  case  relating  a  load  on  a  flat  punch  to  the 

radius  (R)  of  the  base  of  a  cone  crack  propagated  from  the  contact  leading  to 


where  K  is  the  critical  ,-;l  ress  intonsiL''  factor  ai.il  B  is  a  dimension 
less  constant  uniquely  determined  b  the  Poisson's  ratio  (  )  so  that, 
for  a  particular  material,  P  ^  R  ^  .  ilon*  crack  diameters  five  times 
greater  than  the  punch  diameters  were  necessary  to  assure  far  field 
conditions. 

v 

Analogous  P  ■■  c  relations  are  frequently  applied  in  the  case 

in  whicli  the  contact  is  elastic-plastic,  for  example  in  contacts  involving 

(  >  q  i  o ) 

sharp  indenters  ‘  ’  ’  .  Other  types  of  cracks  including  median, 

radial  and  lateral  cracks  are  formed  reflecting  t  he  different  stress 

(2) 

distributions  in  this  case.  (.awn  and  Fuller  have  shown  licit  wedging 
forces  and  friction  can  modify  the  crack  lengths,  with  the  wedging  forces 
depending  on  the  indenter  angle.  Also,  residual  stresses  at  contacts 
can  cause  crack  growth,  especially  near  the  surface  where  the  sign  of 
the  out-of-piane  hoop  stresses  may  be  reversed  from  compressive  to 
tensile  on  unloading^^  . 

i 

In  the  case  of  line  contact  loading,  P  n  l.c  where  L  is  the  length 
of  contact.  Tt  is  clear  that  in  other  cases  in  which  applied  loads  or 
residual  stresses  act  over  an  extended  region,  a  similar  relation  must 
control  crack  propagation  near  the  contact. 


the  basis  that  when  shear  failure  occurs  under  the  contact,  the  out-of¬ 
plane  hoop  tensile  stresses  are  relaxed  to  zero  in  the  plastic  zone.  In 
this  case  the  function  representing  the  radial  dependence  of  the  prior 
out-of  plane  hoop  stress  is  substituted  into  a  standard  fracture  mechanics 
equation  and  integrated  over  the  region  from  the  irreversibly  deformed 
zone  (ZQ)  to  the  tip  of  a  median  crack  of  depth  (c  )  leading  to 
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assuming  c  >>  Z  .  For  v  =  0.25 
m  o 


(2) 


P  =  49.5  K  Z  c2 
I o  n 


(3) 
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Furthermore,  based  on  the  theory  o'"  hardness  testing  “  P  so  t*iat , 

for  a  particular  material,  one  might  expect  to  observe  P  "  c^.  There¬ 
fore,  theoretically  at  least,  the  contact  characteristics  have  a  definite 

effect  on  the  load  dependence  of  the  crack  length  in  this  case  leading 

3/o 

to  P  «  c  rather  than  P  «■  c 

If  one  assumes  that  the  irreversibly  deformed  zones  are  geometrically 
self-similar  for  various  loads,  which  is  reasonable  in  many  cases,  other 
contact  dimensions  such  as  the  indentation  radius  (r^)  or  diameter  (d^) 
can  be  substituted  for  Zq  in  equation  (3).  For  spherical  indenters,  the 
assumption  of  self-similarity  is  a  reasonable  approximation  at  high  loads 
where  there  is  substantial  penetration  of  the  surface.  It  should  be 
noted  that  a  wide  range  of  mechanisms  including  line  contact  loading, 

h 

wedging  and  residur.i  stresses  can  also  lead  to  P  <*  d^  c  relations, 
especially  in  the  near  field.  Therefore,  experimental  agreement  with 
p  «  does  not  necessarily  identify  the  mechanism  controlling  crack 


propagation. 


II.  PROCEDURES 


The  HP  ZnS  was  purchased  as  a  cylindrical  plate  3  in.  diameter  x 
6  mm  thick.  This  plate  was  indented  on  the  large  surfaces  which  are  per¬ 
pendicular  to  the  hot  pressing  direction,  by  static  loading  using  spherical 
indenters.  The  indenters  included  relatively  large  compliant  indenters 
of  glass+  (3  mm  diameter)  and  smaller,  more  rigid  indenters  of  tungsten 
carbide^  with  diameters  ranging  from  0.36  to  1.59  mm.  The  loads  which 
ranged  from  9-890N  were  applied  and  measured  using  a  standard  testing 
machine ^  and  a  crosshead  speed  of  2  ym*s-1.  The  HP  ZnS  plate  was  placed 
on  the  load  cell  and  the  spheres  were  fastened  to  the  crosshead. 

The  damage  was  characterized  by  optical  and  scanning  electron  micro¬ 
scopy.  The  indentation  radius  and  the  length  of  the  longest  radial  crack 
at  each  indentation  were  measured  along  the  specimen  surface  using  the  grid 
in  the  eyepiece  of  the  optical  microscope.  Polarized  light  was  used  to 
improve  the  visibility  of  the  tips  of  the  cracks  during  the  crack  length 
measurements.  Lateral  cracking  and  chipping  around  the  indentations  were 
also  characterized. 

The  results  were  plotted  as  suggested  by  the  relations  P  «  c^2  and 
?  *  di  cr  where  cr  is  the  radial  crack  length.  Coefficients  of  deter¬ 
mination  were  calculated  and  used  to  evaluate  these  relations.  The 
results  are  discussed  in  terms  of  the  similarities  and  differences  in  the 
stress  distributions  beneath  blunt  and  sharp  indenters. 

*  IRTRAN  2,  Eastman  Kodak  Co.,  Rochester,  N.Y.  14650 
+  Glass  spheres,  No.  3000,  Walter  Stern,  Inc.  Port  Washington,  N.Y. 
t  Tungsten  carbide  spheres,  Ultraspherics,  Inc.,  Saulte  Ste.  Marie,  MI.  49783 
$  Instron  Corp.,  Canton,  Mass.  02021 


III.  RESULTS  AND  DISCUSSION 


Relationship  of  load  and  radial  crack  length 

In  the  range  of  loads  used  in  this  investigation,  the  spherical  in- 
denters  induce  permanent  indentations  in  the  HP  ZnS.  Radial  cracks, 
lateral  cracks  and  chipping  were  observed  at  the  indentations.  The 
lengths  of  the  radial  cracks  and  lateral  cracks  increased  with  increasing 
load.  Evidence  presented  previously  includes  descriptions  of  these 
features  including  the  radial  cracks^^’^.  Among  other  characteristics, 
these  cracks  originate  at  the  edge  of  the  indentation  where  they  form 
as  shallow  (Palir.quist)  cracks;  that  is,  as  radial  cracks  that  do  not 
penetrate  very  far  below  the  surface.  As  the  indentation  spreads,  it 
covers  the  original  segments  of  these  cracks.  The  fact  that  these  cracks 
are  observed  within  the  indentation  which  is  in  compression  indicates 
that  these  cracks  grew  during  loading  when  each  segment  was  not  yet  under 
the  indenter.  At  higher  loads  the  cracks  extended  to  greater  depths  so 
that  they  approximated  quarter  penny-shaped  cracks  surrounding  the 
irreversibly  deformed  zone.  Unlike  cracks  formed  under  Vickers  indenters 
where  the  edges  control  crack  locations,  each  individual  crack  does  not 
usually  extend  on  both  sides  of  the  indentation. 

As  shown  previously the  number  of  cracks  formed  increases 
strongly  with  increasing  load,  approaching  20  at  high  loads.  A  wide 
range  of  crack  lengths  is  observed  at  each  indentation. 

Figure  1  shows  that,  for  the  six  glass  and  tungsten  carbide  indent ers, 
the  average  indentation  pressure  remains  approximately  constant  with  in¬ 
creasing  load  indicating  the  relative  absence  of  work  hardening  in  the 
HP  ZnS.  This  result  is  consistant  with  recent  observations  oi  (.hiang, 
Marshall  and  Evans for  ZnS  of  unspecified  type  and  conflicts  with 
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Figure  1*  Indentation  radius  vs.  P  for  spherical  indenters  of  various 
radii  (HP  ZnS) . 
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(4) 

earlier  observations  of  Evans  and  Wilshaw  for  CVD  ZnS.  Work  hardening 

would  be  indicated  by  higher  loads  than  predicted  by  the  linear 

relation,  at  the  larger  indentations.  The  present  average  indentation 

_2 

pressure  is  approximately  1.9  GNm  which  is  consistant  with  the  results 

(16) 


of  Chiang,  Marshall  and  Evans 


for  fully  developed  indentations. 
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The  data  are  plotted  in  terms  of  the  relations  P  ^  c  ' z  and 

r 

U 

P  •*  d^  c^_  where  cr  is  the  length  of  the  longest  radial  crack  in  Figures 
2  and  3,  for  the  two  largest  indenters  which  were  3  mm  diameter  glass 


spheres  and  1.59  mm  diameter  tungsten  carbide  spheres.  The  P  "  c 


V, 


results  using  the  3  mm  diameter  glass  spheres  (open  circles  in  Figure  2) 
yield  a  curved  plot  rather  than  a  straight  line  showing  that  the  response 
of  the  material  is  non-linear  in  terms  of  the  P  ,r  c  relation.  This 


r 

result  indicates  that  another  factor  must  be  taken  into  account. 

The  P  «  d^  c  results  are  represented  by  x  in  Figure  2  yield¬ 
ing  a  straight  line.  This  result  suggests  that  this  relation  is  a  good 
model  for  the  variation  of  load  with  radial  crack  length  under  these 
conditions . 

The  solid  line  through  the  data  points  represents  crack  lengths  cal¬ 
culated  using  equation  (3).  To  do  this  calculation  it  was  first  necessary 
to  calculate  Z  .  This  was  done  by  the  method  indicated  by  Conway  and 
Kirchner^  for  glass  using 


where  4  is  a  proportionality  factor  that  varies  from  one  material  to 
another.  ZQ  =  6a  where  a  is  the  indentation  radius  and  6  is  a  dimension¬ 
less  geometrical  factor  determined  by  the  plastic  zone.  For  circular 
indentations  and  constant  contact  pressure,  substituting  for  (a)  yields 
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Figure  2.  Comparison  of  data  based  on  P  ■  ‘  and  P  •  d.c  relations 

for  3  mm  diameter  glass  spheres  in  HP  ZnS  (c  measured  from 
center  of  indentation).  r 
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Calculated  crack  lengths  assuming  wedging, 
crack  lengths  on  c3/2  scale. 

Calculated  Z0cr,/2  based  on  elastic -plastic 
theory,  equa.  3,  Z0cr,/2  plotted  on 
d;  Cr,/2  scale. 
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Figure  3.  Comparison  of  data  based  on  P  ■■  c  /  and  P  01  d^c  relations 

for  1.59  mm  diameter  tungsten  carbide  spheres  on  HP  ZnS 
(c  measured  from  center  of  indentation). 
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(2irH) 


where  H  is  the  hardness.  Assuming  approximate  self-similarity  and  $ 
constant  for  various  materials  indented  by  spherical  indenters,  C  for 
ZnS  was  estimated  by  multiplying  the  earlier  value  determined  for  glass 
which  was  15  x  10  ^  N  2m, ^  by  the  square  root  of  the  ratio  of  the  hard¬ 
ness  of  glass  (5.6  GPa)  to  the  hardness  of  ZnS  (1.9  GPa)  yielding 
— — W 

25.8  x  10  N  m.  The  Z  values  were  then  calculated  and  substituted 

o 

in  equation  (3)  together  with  KI(,  *  0.44  MPa  m  2  to  calculate  the 

curve  represented  by  the  solid  line.  The  agreement  between  the  data 
points  and  the  calculated  curve  is  remarkable  considering  the  crudeness 
of  the  assumptions. 

A  second  theoretical  curve  indicated  by  the  dashed  line  was  calcu- 

(2) 

lated  using  the  method  of  Lawn  and  Fuller  based  on  the  assumption  that 
wedging  controlled  the  extent  of  crack  growth.  Because  the  wedging 
angle  f  of  a  sphere  is  not  a  fixed  value,  it  was  assumed  in  each  case 
that  the  wedging  angle  was  the  angle  between  the  tangent  to  the  sphere  at 
the  edge  of  the  contact  and  the  perpendicular  to  the  surface  for  each 

indentation.  This  assumption  tends  to  maximize  the  calculated  crack 

3  / 

length.  The  crack  lengths  are  plotted  on  the  c  scale  and  by  comparison 
with  the  experimental  data  (open  circles)  show  that  the  calculated  crack 
lengths  are  much  smaller  than  the  experimental  values.  This  result  indicates 
that  the  experimental  crack  depths  are  not  controlled  by  wedging  in  this 
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Similar  results  are  given  in  Figure  3  for  the  1.59  mm  diameter 

tungsten  carbide  indenter.  Again,  the  experimental  data  plotted  on  the 

3/  % 

cr  2  scale  are  nonlinear.  The  experimental  data  plotted  on  the  d^  c^ 

scale  are  essentially  a  straight  line  indicating  good  agreement  with 

p  cr2'  calculated  crack  lengths  based  on  equation  (3)  yielded  a 

straight  line  with  the  calculated  crack  lengths  slightly  smaller  than 

the  experimental  values,  and  the  wedging  predictions  yielded  smaller 

crack  lengths  than  the  experimental  values,  as  before. 

3/ 

Comparing  the  P  <*  2  curves  in  Figures  2  and  3  reveals  that  in 

each  case  there  is  a  non-linear  region  near  the  origin  but  that,  as  the 

sphere  radius  decreases,  the  region  of  greatest  non-linearity  moved 

closer  to  the  origin.  This  trend  continued  at  still  smaller  sphere 

radii  so  that  at  very  small  sphere  radii  these  curves  form  reasonably 

straight  lines,  despite  the  fact  that  they  tend  to  intersect  the  load 

axis  at  positive  loads.  Therefore,  as  the  spheres  tend  toward  sharp 

3/ 

indenters,  the  P  <*  c  7  relation  becomes  a  somewhat  better  model  than 

r 

it  is  in  the  case  of  indenters  with  larger  radii. 

Based  on  the  above  comparisons,  the  experimental  results  obtained 
using  indenters  with  six  different  diameters  are  given  in  terms  of 

h 

P  <*  d^  c^  in  Figure  4.  The  results  show  approximately  linear  variation 
of  d^  c^  with  increasing  load  and  slightly  decreasing  slopes  (increasing 
radial  crack  length)  with  decreasing  sphere  diameter.  Apparently,  as  the 
sphere  diameter  decreases  there  is  an  increasing  contribution  to  .Tack 
growth  due  to  another  mechanism  such  as  wedging  or  residual  stresses. 

In  some  cases  in  which  there  were  data  points  indicating  less  crack 
growth  than  expected  based  on  the  general  trend  of  the  data  for  a 
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particular  sphere  radius,  extensive  chipping  was  observed  at  the 
indentations.  In  these  cases,  the  residual  stresses  may  have  been 
relieved  to  some  degree  by  chipping  so  that  they  contributed  less  to 
crack  growth.  In  relation  to  the  possibility  of  an  added  wedging 
contribution  to  crack  growth,  it  should  be  pointed  out  that  the  ratio 
of  the  maximum  contact  radius  to  the  sphere  radius  increased  from 
0.26  for  the  3mm  spheres  to  0.77  for  the  .36mm  spheres  so  that,  in 
general,  the  contact  angle  varies  with  sphere  radius  in  the  direction 
expected  to  result  in  an  increased  wedging  effect. 

Statistical  comparison  of  the  P  <*  c  2  and  P  a  d ^  c^2  relations 

The  suitability  of  the  above  relations  was  evaluated  by  comparing 
the  coefficients  of  determination.  When  these  measures  were  originally 
calculated,  the  scatter  introduced  by  errors  in  the  measurement  of  d^  was 
observed  to  effect  the  results.  The  d^  values  are  subject  to  error  because 
the  indentations  are  small  and  not  well  defined,  in  part  because  of 
chipping.  Therefore,  the  coefficients  of  determination  were  recalculated 
using  smoothed  d^  values  from  Figure  1.  These  results  are  given  in 
Table  I.  The  coefficients  of  determination  show  that  the  data  plotted 

h 

according  to  the  P  «  d^  c^,  relation  are  essentially  linear  and  there 
is  little  scatter  about  the  fitted  line.  The  calculated  intercepts 

of  the  six  lines  on  the  load  axis  are  small,  ranging  from  +34  to  -5.4  N. 

3/  , 

The  coefficients  of  determination  for  P  “  c^  *  show  poorer  correlations. 
Also,  the  calculated  intercepts  are  larger,  consistantlv  positive  and 

more  scattered  ranging  from  10.4  to  178.6  N.  These  larger  intercepts 

are  not  unexpected  based  on  the  strong  (almost  parabolic)  curvature 

evident  in  Figures  2  and  3.  Considering  all  of  this  evidence,  it  is 

clear  that  the  data  fits  the  P  a  d^  cr2  relation  better  than  the 

P  «  c  ^  2  relation, 
r 


TABLE  I 


Indenter 

Description 


Coefficients  of  determination  of  data 


3/  u 

in  terms  of  the  P  <*  c  2  and  P  ®  d^  c  models 

r  i  r 


Coefficient  of  Determination 
- f  rf- 


P  «  d.  c  model  P  c 
i  r  r 


3  mm  glass 

1.59  mm  tungsten  carbide 


0.995 

0.997 


0.94 

0.97 


61 


Effect  of  the  stress  distribution  on  crack  propagation 

The  fact  that  radial  cracks  have  been  observed  to  form  on  loading 
(4) 

in  zinc  sulfide  means  that,  at  the  surface,  the  out-of-plane  hoop 

stresses  in  the  elastic  stress  field  are  tensile  under  these  elastic- 

plastic  conditions  rather  than  compressive  as  indicated  by  the  theoretical 

elastic  stress  distribution  for  blunt  indenters  (Figure  5).  Evans  and 
(4) 

Wilshaw  have  shown  numerically  that,  for  the  case  of  a  sphere  pressed 

into  the  surface  under  elastic-plastic  conditions,  to  the  point  that 

the  indentation  radius  equalled  the  sphere  radius,  the  out-of-plane 

hoop  stresses  are  tensile.  Given  the  fact  that  these  stresses  are  the 

result  of  a  central  load  accompanied  by  an  elastic-plastic  disturbance, 

it  is  reasonable  to  assume  that  this  stress  field  varies  inversely  as 

the  square  of  the  radius  from  the  center  of  the.  contact  as  in  the  case 

of  the  out-of-plane  hoop  stresses  directly  under  the  indenter ^ . 

Furthermore,  because  the  magnitudes  of  these  stresses  are.  fixed  by  the 

yield  stress  at  the  elastic-plastic  boundary  and  the  elastic  stress 

(4) 

fields  scale  with  the  radius  of  the  yield  zone  ,  it  is  reasonable  to 
assume  that  the  stress  distributions  along  the  surface  and  under  the 
contact  are  approximately  equal  in  magnitude  and  vary  similarly  except 
to  the  extent  that  the  stresses  along  the  surface  are  affected  by  free 
surface  effects  and  friction.  Therefore,  it  is  reasonable  to  apply 
equation  (3)  to  the  present  radial  cracks  in  addition  to  the  median 
cracks  for  which  it  was  derived. 

As  shown  in  Figure  5,  near  the  surface  under  a  blunt  indenter 
subject  only  to  elastic  deformation,  the  stress  is  compressive  as  one 
would  expect  because  ring  cracks  usually  form  instead  of  radial 
cracks.  The  sign  of  the  stress  changes  to  tensile  between  one 
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and  two  contact  radii  below  the  surface.  At  greater  depths  the  stress 

f 

passes  through  a  peak  and  then  decreases  slightly  more  than  linearly  until 
it  merges  with  the  inverse  square  distribution  in  the  far  field.  On  the 
other  hand  the  stress  distribution  under  the  sharp  indenter  decreases  in¬ 
versely  with  the  square  of  the  normalized  depth.  At  first  glance  it 
might  seem  that  because  of  the  large  difference  in  these  stress  distribu¬ 
tions  the  shift  from  one  of  these  distributions  toward  the  other  might  lead 
to  large  differences  in  crack  propagation.  However,  in  ZnS  there  is 
substantial  irreversible  deformation  under  all  of  the  spherical  indenters. 
Zq  ranges  approximately  from  one  to  two  contact  radii  in  depth.  If  we 
assume  as  we  did  previously  that  the  out-of-plane  hoop  stresses  are 
relaxed  by  the  irreversible  deformation,  it  is  clear  that  the  stress 
intensity  factors  will  not  be  very  different  because  the  stress  dis¬ 
tributions  are  much  less  dissimilar  at  greater  depths. 

The  nature  of  the  near  field-far  field  problem  cno  he  UUc  -  under¬ 
stood  with  the  aid  of  Figures  5  and  6.  Figure  A  is  a  plot  of  the  load 
versus  the  normalized  crack  length  (c  / r^).  The  data  are  plotted  in  this 
way  because  the  boundary  between  the  near  field  and  the  far  field  is 

usually  expressed  in  terms  of  c^/r^.  It  has  been  common  in  the  case  of 
sharp  indenters  to  consider  that  the  transition  from  the  near  field  to 
far  field  behavior  occurs  at  c^/r^  ~  In  Figure  6  this  transition 

should  be  evidenced  by  a  smaller  load  dependence  of  the  crack  length  in 
the  far  field. 

Beginning  with  Figure  5,  one  can  see  that  at  z/r^  >  6  the  stress 
distributions  due  to  the  two  types  of  indenters  become  identical.  However, 
it  would  be  a  mistake  to  conclude  that  the  resulting  crack  lengths  become 
the  same  at  this  point  because  the  stress  intensity  factor  at  the  crack 


Approx,  contours  of  equal  crack  length. 
3  mm  glass 

1.59  mm  tungsten  carbide 
0.81 


Load  vs.  r  /r,  for  spherical  indonters  or  various  diameter 
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tip  depends  on  these  stresses  integrated  from  Z^  to  c.  Therefore,  the 
effect  of  differences  in  contact  characteristics  can  extend  even  farther 
into  the  specimen.  Turning  to  Figure  6,  it  is  clear  that  there  is  r.o 
change  in  slope  or  discontinuity  in  the  data  obtained  using  the  3  ram  glass 
spheres.  The  data  from  the  smaller  tungsten  carbide  indenters  are  some¬ 
what  scattered  so  that  the  results  are  less  clear  but  there  appears  to 
be  a  gradual  change  in  slope  at  c^/r^  -  6  which  may  indicate  a  transition 
from  near  field  to  far  field  behavior.  This  result  is  consistent  with 
Roesler's  observations  for  cone  cracks  and  with  expectations  based  on 
the  stress  distribution  in  Figure  5. 
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IV.  CONCLUSIONS 


1.  For  HP  ZnS  indented  by  spheres  of  various  radii,  the  radial  crack 
length  data  plotted  based  on  P  «  d^,  Cf  2  yielded  essentially  linear 
plots.  Theoretical  curves  based  on  equation  (3)  yielded  good  estimates 
of  the  radial  crack  lengths. 

3/. 

2.  Data  obtained  using  the  larger  indenters  and  plotted  based  on  P  “  c^  x 
yielded  plots  with  substantial  nonlinearity  showing  that  this  is  not  a 

good  model  of  the  response  of  the  material  in  these  cases. 

w 

3.  Small  decreases  in  the  slopes  of  the  P  vs.  d^  c^  curves  with 
decreasing  indenter  radius  may  occur  as  a  result  of  increasing  contri¬ 
butions  of  wedging  or  residual  stresses. 

4.  The  experimental  evidence  suggests  that  for  the  largest  diameter  spheres 
the  observed  cracks  propagated  in  the  near  field.  For  smaller  diameter 
spheres  both  near  field  and  far  field  propagation  may  be  represented. 
However,  the  P  «  d^  c^  relation  fits  the  data  very  well  in  both  cases. 
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ABSTRACT 


The  extent  of  elastic  recovery  at  spherical  indentations  in  MgF2 
indented  by  static  and  impact  loads  was  analysed.  The  fraction  of  the 
indentation  depth  that  was  recovered  decreased  with  increasing  static 
load.  Similar  results  were  observed  for  impact  loads.  The  kinetic 
energy  losses  increased  with  a  velocity  exponent  of  2.61,  close  to  the 
value  2.6  predicted  by  Hunter  for  elastic  impacts,  despite  substantial 
variations  in  the  contributions  of  various  loss  mechanisms  over  the 
range  of  velocities  for  which  measurements  were  made.  Apparently  the 
losses  are  partitioned  among  the  various  loss  mechanisms  that  derive 
their  energy  from  the  elastic  field. 


I.  INTRODUCTION 


Elastic  recovery  at  elastic-plastic  contacts  is  responsible  for 

residual  stresses  at  indentations  and  controls  the  rebounding  velocities 

(1-4) 

of  impacting  particles  .  The  residual  stresses  are  an  important 
factor  in  the  propagation  of  radial  and  lateral  cracks  and  in  determin¬ 
ing  the  strength  and  erosion  rate.  The  rebounding  velocities  are  sensi¬ 
tive  indicators  of  energy  absorption  during  impact.  For  these  reasons, 
understanding  of  elastic  recovery  is  essential  to  understanding  of 

mechanisms  causing  damage  at  contacts. 

(4) 

Lawn  and  Howes  recently  extended  an  earlier  analysis  of  elastic 
recovery  at  conical  indentations  by  Stillwell  and  Tabor introducing 
the  hardness/elastic  modulus  (H/E)  parameter  for  convenience  in  materials 
evaluation  and  incorporating  a  residual  load  term  to  account  for  the 
residual  stresses  necessarily  induced  as  a  result  of  the  post  indentation 
configuration  that  departs  significantly  from  the  initial  flat  surface. 
However,  in  many  important  cases,  contact  damage  is  induced  by  blunt 
rather  than  sharp  indenters.  Therefore,  in  the  present  investigation 
the  Lawn  and  Howes  approach  has  been  used  in  an  attempt  to  analyse 
elastic  recovery  at  spherical  indentations.  This  problem  is  complicated 
by  the  fact  that,  for  spheres,  the  contact  angle  varies  over  the  surface 
of  the  contact  rather  than  remaining  constant  as  in  the  case  of  cones  or 
pyramids.  This  necessitates  the  use  of  assumptions  the  significance  of 
which  must  be  evaluated  by  comparing  predictions  with  experimental  results. 
The  present  problem  is  further  complicated  by  the  fact  that,  for  spherical 
indenters,  the  initial  response  of  the  Indented  material  is  elastic  in 
the  low  load  range  instead  of  elastic-plastic  as  in  the  case  of  sharp 
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indenters.  Therefore,  a  model  that  accounts  for  only  the  elastic- 
plastic  response  can  be  expected  to  yield  good  approximations  only  at 
higher  loads  where  the  initial  elastic  response  is  only  a  small  frac¬ 
tion  of  the  total.  However,  at  these  higher  loads,  radial  and  lateral 
crack  formation,  crushing,  chipping  or  pore  formation  under  the  inden¬ 
tation  may  disrupt  the  simple  indentation  geometry  observed  at  low 
loads . 

Despite  the  fact  that  Lawn  and  Howes  restricted  their  experiments 
to  loads  less  than  the  threshold  for  crack  formation,  it  is  clearly  of 
interest  to  determine  the  applicability  of  their  approach  to  greater 
loads  such  as  those  important  in  investigations  of  localized  impact 
damage  and  erosion.  This  is  one  objective  of  the  present  investigation. 


II.  ANALYSIS 

When  a  rigid  spherical  indenter  is  pressed  into  a  flat  surface  of 
an  elastic-plastic  material,  the  initial  deformation  is  elastic  but  at 
some  load  plastic  deformation  initiates  within  the  plate,  spreads  with 
increasing  load  until  it  reaches  the  sphere-plate  interface  and  then 
along  that  interface  until  it  reaches  the  edge  of  the  contact  or  beyond^. 
The  loading  tends  to  depress  the  edge  of  the  contact  below  the  original 
surface  of  the  plate  but,  depending  on  friction  and  other  factors, 
plastic  flow  toward  the  surface  may  form  a  ring-like  hump  surrounding 
the  contact.  Although  the  profile  of  this  hump  can  be  measured  aitur 
unloading,  little  is  known  about  the  time  dependence  oi  this  process 
during  static  and  impact  loading.  Therefore,  it  is  assumed  in  the  present- 
case  that  these  effects  offset  each  other  so  that,  as  ^hown  in  Figure  la, 


a  Sphere  indenting  flat  surface 


b.  Surface  after  unloading  showing  elevation 
caused  mainly  by  lateral  cracking 


Figure  1  Indentations  during  loading  and  after  unloading 
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the  depth  of  the  indentation  Z  can  be  measured  from  the  original 
flat  durface. 

Loading 

Assuming  as  Lawn  and  Howes  did  that  the  loading  half  cycle  is 

mainly  an  elastic-plastic  process,  the  load  P  is 

P  =  p  -irr .  2  (1) 

o  1 

where  p^  is  the  average  pressure  on  the  projected  area  of  the  indenta¬ 
tion  and  r^  is  the  radius  of  the  indentation.  However,  r?'  =  2RZ  -  7?  -■ 
2RZ  for  shallow  indentations  where  R  is  the  sphere  radius  so  that 

P  =  2 ir  H  R  Z  (2) 

where  the  hardness  (H)  is  taken  equal  to  p^. 

Unloading 

During  unloading  the  material  in  the  plate  pushes  against  the 

indenter  so  that  elastic  energy  is  recovered  from  the  plate.  However, 

as  unloading  progresses,  residual  stresses  are  induced  in  the  region 

near  the  indenter  as  the  elastically  deformed  material  relaxes  against 

the  irreversibly  deformed  material  under  the  indentation.  Therefore, 

only  part  of  the  stored  elastic  energy  is  recoverable  during  unloading. 

During  subsequent  loading  cycles,  the  load  path  is  retraced  with  only 

*  (4) 

small  deviations  as  long  as  the  maximum  load  P  is  not  exceeded  .  Tn 
the  absence  of  cracking  and  other  complications,  unloading  can  be  con¬ 
sidered  to  be  mainly  an  elastic  process. 

As  indicated  by  Johnson^,  the  stresses  in  the  surface  of  an 
elastic,  semi- inf inite  solid  indented  by  a  rigid  sphere  are  (for  small 
strains) 
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_4_  E(yR> 

°  3  TT  (1-V2) 

where  v  is  poisson’s  ratio.  Therefore  the  load  is 

p  _  E  R1^2  Z3^2 

3  (1-v2) 

in  the  absence  of  residual  stresses.  However,  at  each  point  during 
recovery 

P  =  PE  -  Pr 

where  P  and  P  are  the  elastic  and  residual  loads  respectively. 

Cj  r 


(3) 


(4) 


(5) 


Therefore, 


p  =  E  R1/2  (z3/2_  z  3/2) 

3  (1-v2)  ~r 


(6) 


At  the  end  of  the  loading  cycle  and  the  beginning  of  the  unloading  cycle 
the  indenter  loads  (P*)  and  depths  (Z*)  are  necessarily  equal  so  that 


equating  (2)  and  (6), 


Z*"^-  Z  ^2 

r 


Z* 


/6  71  (1-v2)  R1/2\  H 

'  8V2  '  k 


and  assuming  v  =  0.3 


7*3/2.  7  3/2 

- - - *  1.52  R1/2  £ 

Z*  E 


(7) 


(8) 


z*3/2_z  3/2 

so  that  as  one  might  expect  when  H  is  large  relative  to  E,  J _ r_ 


is  large  and  Z^  is  small  (shallow  indentation).  If 


Z* 


Z*  and  Z^.  can  be  measured  or  estimated  it  should  be  possible  to  estimate 


the  dynamic  hardness  during  impact  using  equations  (7)  or  (8). 


Kinetic  energy 

In  the  impact  case,  the  kinetic  energy  is  reduced  to  zero  during 
the  loading  half  cycle  so  that  the  work  done  on  the  plate  (W^)  is 


r 


» 


i 


W.  =  /  P(Z)  dZ 

l 

o 

W.  =  ttHRZ* 
i 


Therefore,  comparing  the  kinetic  energy  of  the  sphere  with  the  work 
done  during  impact  is  another  means  to  evaluate  this  approach. 


III.  EXPERIMENTAL  RESULTS  AND  DISCUSSION 

An  MgF^  plate^  was  subjected  to  static  and  impact  loading  using 

+ 

1. 59  mm  diameter  tungsten  carbide  spheres  as  described  in  more  detail 
elsewhere ’ru~ - *■-“*•  A — • —  ~u - . ~-J  k,.  — 1  - - 


The  contact  damage  was  characterized  by  optical  micro¬ 


scopy  and  prof ilometry .  Good  indentations  were  observed  at  all  static 
loads  in  the  load  range  investigated.  Except  at  the  lowest  loads, 
radial  and  lateral  cracks  were  observed.  In  the  impact  specimens,  the 
lateral  crack  lengths  were  greater  than  the  radial  crack  lengths.  In 
many  cases  these  cracks  extended  to  the  surface  causing  extensive 
chipping.  Because  these  chips  include  a  segment  of  the  surface  of  the 
indentation,  extensive  chipping  interfered  greatly  witli  determination 
of  the  indentation  dimensions. 


Static  loading 


Based  on  measurements  of  the  indentation  radii  after  unloading, 

the  average  pressure  on  the  projected  area  of  the  indentation  is 

approximately  constant  for  the  range  of  loads  used  as  shown  in  Figure 

2  and  Table  I.  The  average  of  these  average  pressures  was  7.86  GPa. 

This  value  is  slightly  higher  than  available  hardness  values,  for 
(4) 

example  7.3  GPa  ,  for  this  material.  These  constant  pressures  support 
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Indentation  Depths  at  Maximum  Load  and  After  Recovery 
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the  assumption  of  constant  pressure  as  a  function  of  Z  in  equation  (2), 
at  least  for  this  particular  range  of  loads. 

The  indentation  depth  at  each  maximum  load  (Z*)  was  calculated 
from  the  recovered  indentation  radii  assuming  that  the  indenter  remained 
spherical  and  using  E  =  115  GPA. 

The  recovered  indentation  depths  (Z  )  were  measured  by  profilometry . 
These  results  are  also  presented  in  Table  I  and  show  that  the  recovered 
depth  ranges  from  about  1/3  to  2/3  of  the  depth  at  maximum  load  and  that 
the  fractional  extent  of  recovery  decreases  with  increasing  load.  This 
variation  is  different  from  the  constant  Z^/Z*  predicted  by  Lawn  and 
Howes  for  pyramidal  indentations  in  particular  materials  and  it  appears 
to  arise  from  the  lack  of  a  constant  indenter  angle  in  the  case  of 
spherical  indentations. 

The  theoretical  predictions  and  experimental  data  for  the  fraction 
of  Z*  recovered  are  compared  in  Figure  3  and  show  reasonable  agreement. 

At  low  values  of  contact  radius  the  hardness  is  not  exceeded  at  any  point 
under  the  contact  so  an  indentation  is  not  observed.  However,  as  the 
load  and  contact  radius  increase,  the  maximum  pressure  will,  at  some  point, 
exceed  the  hardness  under  the  center  of  the  indentation.  If  one  assumes 
that  a  permanent  indentation  can  be  formed  at  this  load  it  is  clear  that 
the  radius  of  such  an  indentation  would  be  small  relative  to  the  total 
contact  so  that  the  assumptions  relating  the  load  and  the  average  pressure 
in  equation  (1)  and  substitution  of  the  hardness  for  the  average  pressure 
in  equation  (2)  are  not  justified  at  low  loads.  The  arrows  in  Figure  3 
indicate  the  Z*  values  at  which  the  maximum  pressure  at  the  center  of 
the  contact  and  the  mean  pressure  reach  7.86  GPa.  Because  permanent 
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Depth  at  Max.  Load  (Z*)-;tm 


f  Figure  3  Fraction  recovered  vs  depth  at  maximum  load  for  static 

indentations  in  MgF2  (WC  sphere  1.59  mm  diam.) 
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indentations  may  be  induced  for  2*  >  4  pm,  it  is  clear  that  equation  (8) 
overestimates  the  fraction  recovered  in  the  range  4<Z*<8.6  pm.  How¬ 
ever,  the  experimental  data  may  extrapolate  to  total  recovery  at 
Z*  =  4  pm,  where  plastic  deformation  is  first  expected  to  occur. 

Impact  loading 

Except  at  low  impact  velocities,  extensive  chipping  removed  much  of 
the  indentation  surfaces  so  that  the  specimens  were  not  suitable  for 
measuring  the  recovered  indentation  depth.  Therefore  available  data  are 
limited  to  three  data  points  mainly  at  low  velocities.  By  substituting 
Z*  and  Z^  in  equation  (8),  the  dynamic  hardness  was  estimated.  The 
results  are  given  in  Table  1  and  plotted  in  Figure  2.  The  fraction 
recovered  was  also  calculated  and  plotted  in  Figure  3.  Because  the 
elastic  recovery  is  somewhat  less  than  expected  for  a  material  of  this 
hardness  and  Young's  modulus,  the  dynamic  hardness  values  and  the  frac¬ 
tion  recovered  are  lower  than  might  have  been  expected. 

The  kinetic  energies  of  the  impacts  can  be  estimated  from  the  mass 
and  velocity  of  the  Impacting  sphere  as  well  as  from  the  depth  Z*  of 
the  indentation.  Results  from  both  methods  are  compared  in  Figure  4. 

At  low  impact  velocities  the  kinetic  energies,  calculated  using  equation 
(8),  are  close  to  the  values  calculated  from  the  mass  and  impact 
velocity.  In  this  velocity  range  good  indentations  with  little  chipping 
are  formed.  Therefore,  the  measurements  of  the  indentation  radii  and 
the  estimates  of  Z*  should  be  reliable.  Furthermore,  the  elastic- 
plastic  mechanism  of  indentation  formation  and  energy  absorption  is 
consistent  with  the  assumptions  of  the  analysis.  Therefore,  the  model 
can  be  used  to  estimate  the  impact  energy  at  low  velocities.  At 
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Figure  4  Kinetic  energy  vs  impact  velocity  for  impacts  of  1.59  mm 
diameter  tungsten  carbide  spheres  on  MrF9. 
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higher  velocities  considerable  chipping,  lateral  cracking  and  crushing 
are  observed.  The  chipping  introduces  uncertainty  into  the  measurement 
of  the  indentation  radii  which  in  turn  affects  the  estimates  of  Z*.  It 
appears  that  in  some  cases  in  which  the  sides  of  the  indentation  have 

been  chipped  away,  r^  has  been  underestimated  leading  to  underestimates 

2 

of  Z*.  Because  the  kinetic  energies  depend  on  Z*  ,  the  resulting  errors 
can  be  large. 

In  principle  it  should  be  possible  to  derive  equations  for  the 
rebounding  kinetic  energy  and  the  coefficient  of  restitution  using  the 
Lawn  and  Howes  approach.  However,  in  view  of  the  scatter  in  the  impact 
kinetic  energy  data,  this  approach  seems  unpromising  in  this  case.  How¬ 
ever,  the  measured  impact  and  rebounding  velocities  were  used  to  calcu¬ 
late  the  coefficients  of  restitution  as  shown  in  Figure  5.  The  coeffi¬ 
cient  of  restitution  (e)  decreases  gradually  with  increasing  impact 
velocity  in  contrast  to  the  Lawn  and  Howes  prediction  of  a  constant 
value  for  impacts  of  pyramidal  indenters.  This  difference  can  be 
explained,  at  least  in  part,  by  the  fact  that  at  low  impact  velocities 
spheres  are  expected  to  have  high  coefficients  of  restitution  because 
there  is  little  or  no  plastic  deformation.  At  higher  impact  velocities 
the  fraction  of  the  kinetic  energy  loss  due  to  plastic  deformation 
increases,  thus  decreasing  the  coefficient  of  restitution.  It  is 
expected  that  this  tendency  will  be  moderated  to  some  extent  by  the 
increase  in  dynamic  hardness  with  impact  velocity. 

The  log  kinetic  energy  loss  varies  linearly  with  log  impact  velocity 
as  shown  in  Figure  6.  The  slope  is  2.61,  a  value  very  similar  to  the 
value  2.6  predicted  by  Hunter for  elastic  impacts  of  spheres  on  flat 
plates.  Except  perhaps  at  6.25  m*s  ^  impact  velocity,  cracks  were 
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observed  at  all  of  the  impact  sites.  The  threshold  velocities  for 
chipping  (10  ra-s  *)  and  crushing  (23  m-s  S  are  marked  on  the  curve. 

One  might  have  expected  that  additional  loss  mechanisms  such  as  chipping 
and  crushing  would  cause  increasing  kinetic  energy  loss  over  that 
expected  but  this  was  not  observed.  Apparently,  because  these  mechanisms 
depend  on  the  elastic  stress  field,  the  energy  losses  are  partitioned 
among  the  available  mechanisms.  Earlier  observations  involving 

crushing  in  a  viscoelastic  material  gave  a  different  result  in  which  the 
slope  increased  at  the  crushing  threshold,  probably  because  of  friction 
losses  that  do  not  depend  on  the  elastic  field. 

The  results  in  Figure  6  are  also  significant  in  terms  of  elastic 
recovery.  One  might  have  thought,  for  example,  that  relief  of  the 
residual  stresses  induced  during  elastic  recovery  by  chipping  might 
increase  the  extent  of  elastic  recovery,  thus  increasing  the  coefficient 
of  restitution  ind  decreasing  kinetic  energy  loss  in  those  specimens. 
Instead,  the  results  seem  to  indicate  that  the  details  of  the  recovery 
process  are  not  important. 

IV.  CONCLUSIONS 

The  existence  of  other  damage  mechanisms  in  addition  to  plastic 

deformation  at  the  indentation  is  not  in  itself  a  decisive  reason  for 

(4) 

not  applying  the  Lawn  and  Howes  approach  to  analysis  cf  the  recovery 
at  indentations  formed  by  static  or  impact  loading  by  spheres.  In  fact, 
an  analysis  for  spheres  based  on  their  approach  yielded  predictions  that 
were  confirmed  in  several  respects  by  static  and  impact  loading  experi¬ 
ments  on  MgF2>  The  principal  difficulties  were  practical  problems  like 
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measuring  the  indentation  radius  when  part  or  all  of  the  sides  of  the 
indentation  have  been  removed  by  chipping  and  measuring  the  depth  of 
the  indentation  in  the  presence  of  a  raised  specimen  surface  caused  by 
lateral  cracking. 

In  MgF2 ,  the  contact  pressure  was  approximately  constant  over  the 
range  of  loads  used  in  the  static  loading  experiments.  The  fraction  of 
the  indentation  depth  that  is  recovered  decreased  with  increasing  static 
load  as  predicted  by  the  theory.  This  result  differs  from  the  constant 
value  obtained  previously  by  Lawn  and  Howes  for  static  pyramidal  indenta¬ 
tions  in  a  wide  variety  of  materials. 

In  impact  experiments,  elastic  recovery  is  similar  to  that  in  static 
loading  experiments.  The  coefficient  of  restitution  decreases  with 
increasing  impact  velocity  reflecting  at  least  in  part  the  decrease  in 
the  fraction  of  the  indentation  depth  that  is  recovered,  with  increasing 
impact  velocity.  The  slope  of  the  log  kinetic  energy  loss  vs.  log  impact 
velocity  curve  was  2.61  in  excellent  agreement  with  Hunter's^1 ^  predic¬ 
tion  for  elastic  impacts  where  the  losses  are  mainly  due  to  stress  waves. 
The  presence  of  different  loss  mechanisms  such  as  plastic  deformation, 
chipping  and  crushing  did  not  vary  the  slope  of  this  curve.  Apparently, 
the  energy  losses  for  the  mechanisms  deriving  their  energy  from  the 
elastic  field  are  subject  to  a  partitioning  process  so  that  if  one 
mechanism  requires  more  energy,  other  mechanisms  must  be  limited. 
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ABSTRACT 


Stress  wave  energy  is  the  principal  energy  loss  mechanism  in 
elastic  impacts  on  firmly  supported  targets  and  is  a  significant  loss 
mechanism  in  other  cases.  In  this  paper,  Hunter's  analysis  is  extended 
to  the  case  in  which  the  mechanical  impedance  of  the  impacting  particle 
is  finite.  Expressions  are  derived  for  the  total  stress  wave  energy  and 
the  stress  wave  energies  in  the  target  and  the  impacting  particle.  These 
expressions  are  used  to  explore  the  effect  of  the  impedance  on  the 
efficiency  of  an  impact  and  the  stress  wave  energy.  The  theoretical 
results  are  compared  with  measurements  of  the  coefficient  of  restitution 
for  impacts  on  a  glass. 
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I.  INTRODUCTION 


Some  theories  of  impact  assume  that  the  stress  wave  energy  is 
negligible  This  assumption  enables  one  to  directly  derive  the  velocity 

dependence  of  important  characteristics  such  as  load,  contact  time  and 
depth  of  damage.  However,  if  one  analyzes  the  energy  losses  during  impact 
on  firmly  supported  targets,  the  stress  wave  energy  is  an  essential  factor, 
especially  in  an  elastic  situation  where  it  represents  the  only  signifi¬ 
cant  loss . 

(?) 

The  present  analysis  is  concerned  with  extension  of  Hunter's 

analysis  of  the  stress  wave  energy  in  a  semi-infinite  solid  to  the  case 

of  a  real  impact  in  which  the  mechanical  impedance  of  the  impacting 

particle  is  not  assumed  to  be  Infinite.  This  analysis  is  followed  by  an 

inquiry  into  the  role  of  the  impedance  of  the  target  and  the  particle  on 

the  efficiency  of  an  impact  and  its  relation  to  the  energy  of  the  stress 

waves.  The  theoretical  results  were  compared,  to  the  extent  possible, 

using  experimental  data  for  the  coefficient  of  restitution  for  impacts 

(3) 

at  various  velocities  on  glass  at  various  temperatures  . 
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II.  EXTENSION  OF  HUNTER'S  ANALYSIS 

(2) 

Hunter  derived  an  expression  for  the  stress  wave  energy  stored  in 

a  target  after  impact  by  a  spherical  projectile.  The  derivation  assumed 

* 

the  radiation  impedance  of  the  projectile  to  be  infinite  in  the  same  way 

(4  5) 

as  Miller  and  Pursey  *  assumed  a  perfect  oscillator  in  their  analysis 

of  the  field  and  radiation  impedance  of  mechanical  radiators. 

However,  in  real  impact  situations,  the  radiation  impedance  of  an 

impacting  particle  (i.e.  the  radiator)  is  finite.  If  a  stress  is  applied 

to  the  particle,  a  displacement  will  appear  within  it,  and  stress  waves 

will  propagate.  Thus,  part  of  the  stress  wave  energy,  produced  by 

impact,  is  induced  in  the  particle. 

To  derive  an  expression  for  the  stress  wave  energy  induced  in  the 

particle,  one  can  follow  a  derivation  similar  to  that  of  Miller  and 
(4) 

Pursey  for  a  semi-infinite  solid.  This  approach  is  reasonable.  Indeed, 

Love^  has  shown  the  similarities  in  the  elastic  behavior  of  a  half 

plane  and  a  disk.  Therefore,  one  can  expect  such  similarities  to  come 

into  play  when  comparing  a  sphere  and  a  semi-infinite  solid. 

(4) 

Miller  and  Pursey' s  derivation  describes  the  case  where  a  circular 
disk  of  finite  radius  vibrates  normally  to  the  free  surface  of  a  semi¬ 
infinite  solid.  For  these  conditions,  they  derive  the  expression  for  the 
displacement  of  the  free  surface  of  the  solid.  Then  they  determine  the 
energy  carried  by  the  stress  waves  thus  produced  and  compute  the  partitioning 

*The  radiation  impedance  is  defined  as  the  ratio  of  the  stress  to  the 
mean  displacement  velocity  under  the  radiator  in  the  direction  of  the 
applied  stress^). 


-k. 


of  this  energy  among  the  various  types  of  waves  (compressional,  shear  and 

surface  waves)  Hunter  generalized  the  results  of  Hiller  and  Pursey 

to  cover  the  problem  of  transient  pulses .  Using  the  expression  for  the 

normal  displacement  of  the  free  surface  of  a  semi-infinite  body  subjected 

to  a  uniform  periodic  surface  pressure  pe*Ut  acting  over  a  surface  of 

radius  (a) ,  the  mean  surface  displacement  (u)  over  the  circle  of  contact 

(contact  area  of  impact)  was  obtained,  by  Fourier  synthesis,  for  a  pressure 

(2) 

pulse  of  arbitrary  profile  p(t)  .  This  result  was  a  key  to  determining 

the  stress  wave  energy  produced  in  the  target  during  impact  which  is 

given  by*  2  T  m  -1/5  R  R“6/5 

p  P  3/5  2 

Wt  = - 3 - V0  mp  V  (1) 

pt  C0t 


For  a  spherical  particle,  the  derivation  will  follow  the  same  pattern 


described  in  the  above  paragraph 

.(4) 


As  Miller  and  Pursey  did,  oue  can  write  the  expression  relating 
the  components  of  the  displacement  (u^)  to  the  dilatation  (div  5  =  A)  and 
the  distortion  (rot  U  -  W) ,  as  well  as  the  partial  differential  equations 


characteristic  of  the  dilatation  and  the  distortion.  The  expressions  for 


the  stress  across  the  free  surface  as  functions  of  the  displacement  complete 


$  (7 

the  set  of  equations  necessary  to  solve  the  problem.  Hankel's  transform  ' 


was  used  to  change  the  partial  differential  equations  to  differential 
equations.  In  the  case  of  a  semi-infinite  solid,  this  transformation  is 


See  Appendix  A  for  definitions  of  symbols. 


These  relationships  are  written  in  a  cylindrical  system  of  coordinates 
with  the  origin  at  the  center  of  the  contact  area  and  the  z  axis  extending 
from  this  origin  along  the  diameter  of  the  particle  (Fig.  1)  as  in  the 
case  (z)  of  Reference  4. 


Figure  1.  System  of  coordinates  used  in  the  derivation  (center  is 

at  the  center  of  the  contact  area.  System  (0,r)  parallel 
to  the  target  and  axis  7.  coinciding  with  the  diameter  of 
the  particle  and  with  origin  0). 
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relatively  simple  but  complications  occur  in  the  spherical  case  because 

application  of  Hankel's  transform  to  a  derivative  involves  introduction 

of  boundary  conditions.  These  boundary  conditions  are  such  that  the  surface 

of  the  sphere  is  free  of  stress  outside  the  contact  area.  Under  these 

conditions,  remembering  cylindrical  symmetry,  use  of  Hankel's  transform 

implies  introduction  of  bulk  dilatation  and  bulk  distortion  of  the  sphere 

and  their  derivatives.  In  effect,  in  the  treatment  of  the  stress  waves 

in  a  sphere,  the  partial  differential  equations  are  transformed  into  non- 

homogenous  differential  equations  (see  Appendix  B) .  These  equations 

could  probably  be  solved  but  since  the  analysis  of  the  bulk  deformation 

of  a  particle  is  beyond  the  scope  of  this  particular  investigation,  this 

problem  was  avoided  by  assuming  a  rigid  sphere  (based  on  the  distinction 

that  stress  waves  are  properties  of  the  medium  and  bulk  deformation  is  a 

A 

property  of  the  body  ) .  This  assumption  enables  one  to  equate  to  zero 

the  second  member  of  the  differential  equations  characteristic  of  the 

stress  waves  in  the  sphere  and,  in  consequence,  to  use  the  results  for 

a  semi-infinite  solid  to  describe  these  stress  waves  in  a  sphere.  In 

order  to  determine  the  stress  wave  energy  in  the  particle,  we  can  relate 

it  to  the  amount  of  stress  wave  energy  in  the  target  through  the  notion 

of  intensity.  The  intensity  of  an  elastic  wave  is  defined  as  the  mean 

(9) 

power  transmitted  per  unit  of  surface  .  For  the  target,  this  intensity 


Such  an  assumption  is  comparable  to  the  assumption  of  a  non¬ 
bending  plate. 


will  be  equal  to  the  stress  wave  energy  divided  by  the  time  of  contact 


(tc)  and  the  contact  area  (ira  ; 


W. 


Xt  = 


(Jm“2s_1) 


(2) 


t  7ra 
c 


This  wave  intensity  is  known  to  be  equal  to  the  applied  stress  squared 

(9) 


divided  by  the  characteristic  impedance  of  the  medium  (PtcQt) 


Therefore,  one  can  write  for  the  target 


Xt  =  — 2 
t  ifa 
c 


Vot 


(3) 


Similary,  for  the  particle 


W 


I  = 
P 


t  na 
c 


°pc0p 


f4) 


As  p,  tcand  ira  are  the  same  for  the  target  and  particle,  they  can  he 


eliminated  yielding 


W  =  W  (  — — — ) 

p  c  ppcop 


(5) 


Therefore,  Hunter's  formula  can  be  extended  as 


W 


w  (i  +  -  Pt  --t- 
t  p_c, 


(6) 


P  Op 


in  which  is  given  by  equa.  (1).  This  equation  shows  that  the  impedance 


ratio  between  the  target  and  the  particle  has  a  significant  effect  on  the 
impact  characteristics. 
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III.  EFFECT  OF  THE  IMPEDANCE  RATIO  ON  THE 
PARTITIONING  OF  STRESS  WAVE  ENERGY 


If  one  assumes  the  particle  to  have  infinite  impedance,  one  derives 

(2) 

the  results  found  by  Hunter  ,  namely 

W8  ■  “t  <7) 

If  one  assumes  the  particle  and  the  target  to  have  the  same  character¬ 
istic  impedance,  one  should  find  that  the  total  stress  wave  energy  is  twice 
that  induced  in  the  target.  Kirchner  and  Gruver^*^  obtained  experimental 
results  consistent  with  this  expectation  for  impacts  of  glass  spheres  on 
glass  plates  for  which  the  target  and  particle  have  the  same  characteristic 
impedance.  Similarly,  this  result  implies  a  better  match  between  Hunter's 
theory  and  the  data  of  Tillet^^.  In  that  case  the  experimental  values 
for  the  coefficient  of  restitution  are  lower  than  predicted  by  Hunter's 
theory.  Introduction  of  a  finite  impedance  ratio  implies  a  larger  stress 
wave  energy  loss,  hence  a  lower  coefficient  of  restitution. 

The  distribution  of  the  stress  wave  energy  between  the  particle  and 
the  target  varies  substantially  depending  on  the  properties  of  the  parti¬ 
cular  materials.  Letting  (n)  be  the  ratio  of  the  characteristic  impedance 
of  the  target  to  that  of  the  particle,  equa.  (6)  becomes 


W  =  W_  (1  +  n) 
st 

Wt  1 

W  1  +  n 
s 


W 


W_ 


n 


(8) 

(9) 


and 


1  +  n 


(10) 


I 


Table  I  gives  the  percentages  of  the  total  stress  wave  energy  in  the  target 
and  the  particle  for  the  realistic  range  of  (n) .  These  data  have  practical 
importance  in  some  cases.  For  example,  ceramic  armor  with  a  high  Young's 
modulus  (i.e.  a  high  characteristic  impedance)  is  likely  to  have  a  high 
impedance  compared  to  a  given  projectile  so  that  substantial  stress  wave 
energy  is  induced  in  the  projectile,  increasing  the  likelihood  that  the 
projectile  will  be  destroyed  by  reflected  stress  waves.  Similarly,  in 
the  case  of  IR  windows,  high  characteristic  impedance  is  a  desirable 
property  for  the  same  reason. 


TABLE  I 

PERCENTAGES  OF  THE  TOTAL  STRESS  WAVE  ENERGY  IN  THE 
PARTICLE  AND  TARGET  FOR  VARIOUS  VALUES  OF  n 
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IV.  THE  EFFECT  OF  TARGET  AND  PARTICLE  IMPEDANCE  ON  THE 
AMOUNT  OF  STRESS  WAVE  ENERGY  PRODUCED  DURING  IMPACT 


Before  analyzing  the  expression  for  total  stress  wave  energy 
developed  in  Section  II  (equa.  (6)),  it  is  enlightening  to  analyze  in 
detail  Hunter's  formula  (equa.  (1))  in  terms  of  impedance.  A  simple 
substitution  of  the  value  of  the  target  impedance 


zt  =W’7 =  pt  cot 


in  equation  (1)  yields 


2,«  -1'5  *3'5  s‘6/5  P  2  V  325 
p  to 


b  -p  vo] 


This  expression  cannot  be  reduced  to  a  simple  function  of  the  target 
impedance  *  f(Zt).  Thus,  it  must  be  analyzed  as  a  multivariate 
function  dependant  upon  impedance,  density  and  elastic  modulus  of  both 
target  and  projectile,  or  formally 

“t  -  f(V  V  V  V  “p>  <l: 


To  evaluate  the  behavior  of  such  a  function,  a  study  of  its 

projections  in  various  planes  is  necessary,  a*  First  the  projectile 

9  2 

and  the  incoming  velocity  are  kept  constant  (E^  *  50  10  N/m  ,  = 

3 

2000  kg/m  VQ  ■  80  m/s  R  -  .001m),  thus  insuring  constant  incoming 


energy.  With  this  constraint,  the  density  and  elastic  constant  of  the 
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target  are  independantly  varied  within  a  range  centered  around  the 
projectile  characteristics 

1000  <  <  3000  K  m~3 

t  g 

30 -109  <  E  <  70 -109  N  m~2 

The  variations  of  the  resulting  stress  wave  energy  in  the  target  (W^.) 
are  then  plotted  as  a  function  of  the  density  and  the  elastic  constant 
of  the  target  (Fig.  2a  and  Fig.  2b).  These  graphs  show  a  behavior  consis¬ 
tent  with  simple  wave  mechanics;  an  increase  in  the  density  (i.e.  mass) 
yields  more  stress  wave  energy  in  the  form  of  inertia.  On  the  contrary, 
an  increasing  elastic  constant  decreases  the  potential  for  displacement 
(particle  velocity)  and  thus,  decreases  the  energy  of  the  stress  wave. 

Fig.  2c  shows  the  variation  of  the  stress  wave  energy  in  the  target 
(W^.)  as  a  function  of  the  target  impedance.  The  variations  at  constant 
elasticity  or  constant  density  are  easily  traced  on  this  graph.  This 
plot  also  shows  that  the  spread  of  possible  values  of  the  stress  wave 
energy  varies  with  the  impedance,  as  a  relatively  larger  target  impedance 
yields  less  variability  than  a  smaller  one.  This  effect  is  due  to  the 
relative  importance  of  the  target  density  and  elastic  modulus  in  the 
impedance,  which  favors  the  elastic  modulus  at  high  impedance  and  favors 
the  density  at  low  impedance. 

b.  To  complete  the  analysis  of  the  stress  wave  energy  in  the  target, 
a  similar  evaluation  of  Hunter's  equation  must  be  made  for  a  constant  target. 
However,  to  allow  comparison  with  the  previous  analysis,  one  must  keep 
the  same  incoming  energy.  As  a  consequence,  one  must  vary  the  incoming 
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Figure  2b.  Stress  wave  energy  in  target  vs.  Young's  modulus  of 

target  for  fixed  particle  properties  and  impact  energy. 
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velocity  with  the  varying  particle  density  according  to 


where  V  and  p  are  the  velocity  and  density  of  the  particle  in  the 
o  o 

previous  analysis.  With  this  minor  change,  the  analysis  proceeds  as 
before  (Fig.  3a,  3b  and  3c).  There  the  plot  of  the  stress  wave  energy 
shows  little  variation  with  density  and  a  decreasing  variation  with 
the  elastic  constant.  One  must  remark,  however,  that  the  energy  developed 
in  that  case  is  somewhat  lower  than  in  the  previous  case;  in  particular, 
the  level  of  energy  varies  when  the  characteristics  of  the  target  and 
the  particle  are  interchanged. 

When  the  total  stress  wave  energy  is  considered,  two  cases  must  be 
distinguished  (Fig.  4a  and  4b): 

a.  If  the  particle  characteristics  are  kept  constant  and  the 
target  characteristics  are  taken  as  equal  to  those  of  the  particle, 
then  the  total  stress  wave  energy  is  double  that  in  the  target.  As 
the  target  impedance  decreases,  the  total  stress  wave  energy  becomes 
less  than  double  the  target  stress  wave  energy  and,  conversely,  as  the 
impedance  of  the  target  increases,  the  total  stress  wave  energy  is  more 
than  double  that  of  the  target.  The  overall  effect  is  a  reduction  in 
the  variation  of  spread  of  possible  values  for  the  stress  wave  energy 
for  various  target  impedances  that  were  observed  earlier  (Fig.  2c). 

b.  When  the  target  characteristics  are  kept  constant,  the  effect 
of  the  impedance  ratio  is  reversed:  low  relative  particle  impedance 
yielding  a  ratio  of  total  stress  wave  energy  to  target  stress  wave 
energy  larger  than  2  while  a  high  relative  particle  Impedance  yields 


a  ratio  lower  than  2. 
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Figure  '3a.  Stress  wave  energy  in  target  vs.  particle  density  for 
fixed  target  properties  and  impact  energy. 
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Figure  '3b.  Stress  wave  energy  in  target  vs.  'dung's  modulus  of 

particle  dr  fixed  target  properties  and  inpact  energy. 
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Figure  4b.  Comparison  of  the  stresd  wave  energy  in  the  target  and 
the  total  stress  wave  energy  for  fixed  target  propertio 
and  impact  energy. 
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In  conclusion,  the  analysis  of  the  variations  of  the  total  stress 
wave  energy  due  to  impact  shows  a  very  complex  dependency  upon  target 
and  projectile  characteristics;  in  particular,  these  variations  are  very 
different  when  the  impedance  is  varied  at  constant  elasticity  and  when 
varied  at  constant  density.  In  fact  these  variations  are  opposite. 
Furthermore ,  it  can  he  noted  that  the  asymmetry  of  the  problem  (i.e.  non 
equivalence  of  energy  when  target  and  projectile  characteristics  are 
interchanged)  is  inherent  to  Miller  and  Pursey*'1  ^  derivation.  the 
introduction  of  the  impedance  ratio  tends  to  temper  this  asymmetry  but, 
further  investigations  into  the  mechanics  of  the  contact  should  be  made 
to  permit  a  physical  explanation  of  this  phenomenon. 


V.  COMPARISON  OF  THE  THEORY  WITH  EXPERIMENTS  USING  GLASS 


(3) 

Kirchner  and  Gruver  determined  the  coefficient  of  restitution  for 

impacts  of  glass  spheres,  3mm  in  diameter,  on  firmly  supported  soda  lime 

glass  plates.  They  found  that,  at  room  temperature,  the  coefficient  of 

restitution  remained  above  0.9  over  a  wide  range  of  impact  velocities 

(Figure  5).  As  shown  in  Figure  6,  the  present  theory  which  assumes  only 

stress  wave  energy  losses  yields  similar  results.  The  effect  of  the 

temperature  dependence  of  density,  Young's  modulus  and  Poisson's  ratio, 

(12) 

was  investigated  using  data  of  Duchateau  and  co-workers  .  Assuming 
only  stress  wave  energy  losses,  the  effect  of  the  variations  in  these 
properties  with  temperature  on  the  coefficient  of  restitution  is  small 
(Figure  6).  The  energy  loss  varies  as  the  square  of  the  coefficient  of 
restitution.  The  variation  at  525°C  is  less  than  that  at  440°C  because 
the  Poisson's  ratio  decreases  sharply  around  the  annealing  temperature 
(500  -  550°c)(12). 

These  variations  due  to  the  stress  wave  energy  losses  do  not  account 

for  the  large  decreases  in  coefficient  of  restitution  observed  at  high 

0 

velocities  at  440  and  525  C.  Therefore,  the  explanation  based  on  increased 
(3) 

crushing  remains  the  only  explanation  available.  At  higher  temperatures 
650  and  720°C,  which  are  above  the  softening  point,  the  variation  of  the 
coefficient  of  restitution  with  velocity  is  different  from  that  observed 
at  lower  temperatures.  Lower  coefficients  of  restitution  (higher  energy 
losses)  are  observed  at  low  impact  velocities  but  the  decreases  are  not 
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as  great  at  higher  velocities.  At  these  higher  temperatures,  the  glass 
deforms  viscoelastically  to  form  indentations.  Because  there  is  no 
threshold  velocity  or  yield  point  for  this  viscoelastic  deformation,  the 
effects  of  the  resulting  energy  losses  are  observed  at  all  velocities. 
Therefore,  it  is  reasonable  to  conclude  that  this  change  in  shape  of  the 
coefficient  of  restitution  vs.  impact  velocity  curves  is  the  result  of 
the  introduction  of  the  indentation  energy  loss  mechanism  at  temperatures 
above  the  softening  point. 
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LIST  OF  SYMBOLS 


Radius  of  the  area  of  contact 

Compressional  elastic  constant  of  the  medium 

Shear  elastic  constant  of  the  medium 

Velocity  of  a  longitudinal  wave  in  a  thin  rod  made  of 

the  same  material  as  the  Impacting  particle  (=y  E^/p^) 

Velocity  of  a  longitudinal  wave  in  a  thin  rod  made  of 

the  same  material  as  the  target  (*^E^7p~) 

Young's  modulus  of  the  particle 
Young's  modulus  of  the  target 


2  2  2 

-  a  -  ixr 

-  o2 

— 2-  (eq  (1)) 
P 


-Y2)* 


(eq  (1)) 


gradient  operator 

Hankel's  transform  of  order  m 

Wave  intensity  in  the  particle 

Wave  intensity  in  the  target 

Bessel's  function  of  order  m(meQ) 

uj(pt/ci  i)*5 

“(pt/c44)JS 

Mass  of  the  particle 

Order  of  Bessel's  function 


Index  characteristic  of  impedance  ratio  (Ptcot^ppcOp^ 


APPENDIX  A  (continued) 


Applied  stress 

Set  of  rational  numbers 

Radial  dimension  of  a  cylindrical  system  of  coordinates 
Radius  of  the  particle 
Rotational  operator 

Contact  time  between  target  and  particle 
Vector  displacement 

Radiai  component  of  the  displacement 
Angular  component  of  the  displacement 
Axial  component  of  the  displacement 
Impact  velocity 

Partial  function  for  integration  by  parts 

Stress  wave  energy  in  particle 

Stress  wave  energy  in  target 

Total  stress  wave  energy 

Partial  function  for  integration  by  parts 

Target  impedance 

Particle  impedance 

Axial  dimension  of  a  cylindrical  system  of  coordinates 
Axial  stress  at  the  free  surface  of  the  sphere 
Radial  stress  at  the  free  surface  of  the  sphere 
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APPENDIX  A  (continued) 

2(l-o  )  , 

<ei  <d» 

dilatation  (=  divergence  U  =  div  IJ) 

Hankel's  variable 

Angular  dimension  of  a  cylindrical  system  of  coordinates 
Lame  constant 
Lam£  constant 

Dummy  variable  for  integration 
Density  of  the  particle 
Density  of  the  target 
Poisson's  ratio  of  the  particle 
Poisson's  ratio  of  the  target 

(1.068)5  (1  +  o.)  (4I>6/5  3(0 

t  1-tO^  lb  t 

Angular  frequency  of  a  periodic  phenomenon 
Laplacian  operator 

•) 

=  7  a7  (r  aF>  +  T  72  +  ;T  in  8  cylindrical 
r  do  dz 


§ 


*■ 


system  of  coordinates 
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APPENDIX  B 


EQUIVALENCE  OF  THE  STRESS  WAVE  ENERGY  OF  A 
SEMI-INFINITE  SOLID  AND  SPHERE 


Starting  from  the  equation  of  elasticity 

-> ,  4  ^  a  a2ii 

'  ll  grad  A  -  c44  rot  W  =  pt  ~2 

in  which  U  is  the  vector  displacement 

A  =  div  U  is  the  dilation 

W  =  rot  U  is  the  distortion 


(15) 


cu  =  X  +  2v 


C44  "  P 


linear  combinations  of  the 
Lame  constants  X  and  p 


To  account  for  the  time  dependence,  we  assume  the  dependence  on  time 
to  be  harmonic  so  that  (15)  becomes 


->  St  a  2  -* 

c^  grad  A  =  c^  rot  W  +  pt«  0=0 


(16) 


This  equation  projected  on  the  axes  z  and  r  of  a  cylindrical  system 
gives  two  of  the  equations  necessary  for  the  solution  of  the  problem 
of  stress  waves  in  a  sphere. 

Taking  the  divergence  and  the  curl  of  equation  (16)  yields  two 
additional  equations  respectively  characteristic  of  A  and  W-  Then, 
writing  the  stress  across  the  free  surface  gives  the  last  two  required 
equations.  These  six  equations  are  listed  as  follows: 
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Ti 

IJ 

s 

* 


3A 


3W 


cn  t;  +  c//  T“  +  p*.u  u  =o 
11  3r  44  3z  t  r 


(17) 


„  9A  44  3(rW)  .  2 

11  te  +  T  ~3r  +  pt“  “z  ■  0 


(18) 


2  2 

Cjj  V  A  +  p^oi  A  =  0 


(19) 


17  (r  T7W_))  +  ^7  +  k2  w  =  0  wlch  k2  =  “ Jr; ■ 

riz  \  44 


(20) 


2  - 


2,  2 


44 

P  to2  7r  _  3  ( J.  3(rW).  _  2  32A  32W 

2  3r  r  3r  ’  ~  y  3r3z  „  2 


jg  =  .2  fr  9W.  _  iK^-2)  3  3Av  _  p4  32A 

.2  3r  U  3r '  r  17  (r  17)  ~~2  (21) 

oZ 


'44 


3z 


(22) 


with  p2  -  -1-1. 

C44 


To  solve  these  equations  one  must  eliminate  r.  This  is  done  by  the 
Hankel  transform  defined  as 


gm(5)  =  /Q  f(r)r  Jn(rc)  dr 


(23) 


and  its  inverse 

f(r)  =  /«  8m(?)  5  J  (rO  d£  (24) 

urn  m 

where  J  (re)  is  a  Bessel  function  of  the  mth  order.  Thus,  the  Hankel 
m 

transform  replaces  the  function  f(r)  by  its  transform  gm(C)  independent 
of  r.  One  further  advantage  of  the  Hankel's  transform  is  to  change 
a  differentiation  into  a  multiplication  as  in  the  case  of  Laplace's  or 


t 


I 


Fourier's  transforms,  and  hence,  change  a  partial  differential 
equation  with  two  variables  into  a  differential  equation.  In  fact, 
to  take  the  transform  of  the  derivative  of  f  with  respect  to  r 


r0  r  Vr5)  dr 


one  must  do  an  integration  by  parts  by  writing 


v  =  r  J  (r?) 
m 


dw  =  Hr1  dr 


(26a) 


(26b) 


which  by  differentiation  of  (26a)  ^  and  integration  of  (26b)  gives 


dv  =  c r  J  .(r£)  dr 
m—  i 


(27a) 


w  =  f(r) 


(27b) 


which  combined  with  (26a)  and  (26b)  gives  the  value  of  the  Integral  (25) 


0  Hr1  r  Jm(r°  dr  “  Cf(r)  r  Jm(r°3o  *  C  /0  f(r)  r  Jm-l(r0  dr  (28) 


Vl(f}  =  Cf(r>  r  Jm(r?)]o  -  ?Vl(f) 


As  stated  earlier,  the  Hankel  transform  has  changed  the  derivative 
3f 

'gj'  into  a  multiplication  ?gm_^(f)  but  it  has  also  introduced  the 
boundary  conditions  [f(r)  r  In  the  case  of  a  semi-infinite 


solid  these  boundary  conditions  vanish.  In  the  case  of  a  sphere,  they 
are  a  function  of  the  bulk  displacement  and  bulk  distortion  of  the 
sphere.  These  displacements  are  due  to  the  deformability  of  the  sphere 
which  allows  the  bulk  vibrations  of  the  sphere However,  when 
the  sphere  is  small,  as  in  the  case  >f  particle  impact,  the  frequency 
of  the  bulk  vibrations  is  high  so  that  they  cannot  be  excited  by  the 
relatively  low  frequency  impact  pulse.  therefore,  in  most  cases, 
these  bulk  vibrations  can  be  neglected.  If,  in  addition,  we  assume 
that  no  plastic  deformation  occurs  in  the  sphere  upon  impact,  the 
boundary  conditions  introduced  by  the  Hankel's  transform  vanish  as 
in  the  case  of  a  semi-infinite  solid.  The  derivation  of  the  expression 
for  the  displacement  thus  follows  that  of  Miller  and  Pursey^11^. 
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